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ASYMPTOTIC PROPERTIES OF LINEAR FIELD EQUATIONS IN 

ANTI-DE SITTER SPACE 

GUSTAV HOLZEGEL, JONATHAN LUK, JACQUES SMULEVICI, AND CLAUDE WARNICK 


Abstract. We study the global dynamics of the wave equation, Maxwell’s equation and 
the linearized Bianchi equations on a fixed anti-de Sitter (AdS) background. Provided 
dissipative boundary conditions are imposed on the dynamical fields we prove uniform 
boundedness of the natural energy as well as both degenerate (near the AdS boundary) 
and non-degenerate integrated decay estimates. Remarkably, the non-degenerate esti¬ 
mates “lose a derivative”. We relate this loss to a trapping phenomenon near the AdS 
boundary, which itself originates from the properties of (approximately) gliding rays near 
the boundary. Using the Gaussian beam approximation we prove that non-degenerate 
energy decay without loss of derivatives does not hold. As a consequence of the non¬ 
degenerate integrated decay estimates, we also obtain pointwise-in-time decay estimates 
for the energy. Our paper provides the key estimates for a proof of the non-linear sta¬ 
bility of the anti-de Sitter spacetime under dissipative boundary conditions. Finally, we 
contrast our results with the case of reflecting boundary conditions. 


1. Introduction 

The non-linear Einstein vacuum equations with cosmological constant A, 

(1) Ric[g\=Kg, 

constitute a complicated coupled quasi-linear hyperbolic system of partial differential equa¬ 
tions for a Lorentzian metric g. The past few decades have seen fundamental progress in 
understanding the global dynamics of solutions to ([^. In particular, a satisfactory answer 
- asymptotic stability - has been given for the dynamics of ([^ with A = 0 near Minkowski 
space [I], and the dynamics of (Q with A > 0 near de Sitter space mm, the maximally 
symmetric solution of ([^ with A > 0. Today, the dynamics near black hole solutions of 
([^ for both A = 0 and A > 0 is a subject of intense investigation |3l O [6]. 
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In contrast to the above, the global dynamics of Q with A < 0 near anti-de Sitter space 
(AdS), the maximally symmetric solution of the vacuum Einstein equations with A < 0, 
is mostly unknown. Part of the problem is that in the case of A < 0, the PDE problem 
associated with Q takes the form of an initial boundary value problem.Therefore, even to 
construct local in time solutions, one needs to understand what appropriate (well-posed, 
geometric) boundary conditions are. It also suggests that the global behaviour of solutions 
starting initially close to the anti-de Sitter geometry may depend crucially on the choice 
of these boundary conditions mu- 

Using his conformal field equations, Eriedrich [8j constructed local in time solutions to 
Q with A < 0. (See also [9] for a recent proof for Dirichlet conditions using harmonic 
gauge.) While in general it is quite intricate to isolate the geometric content inherent in 
the boundary conditions imposed (partly due to the large gauge freedom present in the 
problem), there is nevertheless a “conformal” piece of boundary data that does admit a 
physical interpretation. This goes back to the Bianchi equations, 

(2) [Vg]“IUa6cd = 0, 

satisfied by the Weyl-tensor of a metric g satisfying 0- As is well known pQ , the equations 
0 can be used to estimate the curvature components of a dynamical metric g. The 
boundary data required for well-posed evolution of 0 will generally imply a condition on 
the energy-flux of curvature through the timelike boundary. Two “extreme” cases seem 
particularly natural and interesting: The case when this flux vanishes, corresponding to 
reflecting (Dirichlet or Neumann) conditions, and the case when this flux is “as large as 
possible”, corresponding to “optimally dissipative” conditions. While any such boundary 
data for Q will have to be complemented with other data (essentially the choice of a 
boundary defining function and various gauge choices) to estimate the full spacetime metric 
[3, it is nevertheless reasonable, in view of the strong non-linearities appearing in the 
Einstein equations, to conjecture the following loose statement for the global dynamics of 
perturbations of AdS under the above “extreme” cases of boundary conditions: 

Conjecture 1. Anti-de Sitter spacetime is non-linearly unstable for reflecting and asymp¬ 
totically stable for optimally dissipative boundary conditions. 

The instability part of Conjecture was first made in mm in connection with work 
on five-dimensional gravitational solitons. See also [12]. By now there exist many refined 
versions of this part of the conjecture as well as strong heuristic and numerical support in 
its favour [laiiatts]. 

The present paper is the first of a series of papers establishing the stability part of 
Conjecture Here we contribute the first fundamental ingredient, namely robust decay 
estimates for the associated linear problem. 

1.1. Linear field equations on AdS. An important prerequisite for any non-linear sta¬ 
bility result is that the associated linear problem is robustly controlled |16) . In the present 
(dissipative) context, this means that the mechanisms and obstructions for the decay of 
linear waves in the fixed AdS geometry should be understood and decay estimates with 
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constants depending on the initial data available. We accomplish this by giving a complete 
description of the decay properties of three fundamental field equations of mathematical 
physics: 

(W) The conformal wave equation for a scalar function on the AdS manifolcQ 

(3) + 2tt = 0 . 

(M) Maxwell’s equations for a two-form F on the AdS manifold, 

(4) dF = 0 and d F = 0 . 

(B) The Bianchi equations for a Weyl field W (see Definition below) on the AdS 
manifold 


(5) = 

Note that since AdS is conformally flat, equation ([^ is precisely the linearization of the 
full non-linear Bianchi equations ([^ with respect to the anti-de Sitter metric. 

The models (W), (M), and (B) will be accompanied by dissipative boundary conditions. 
In general, to even state the latter, one requires a choice of boundary defining function for 
AdS and a choice of timelike vectorfield (or, alternatively, the choice of an outgoing null- 
vector) at each point of the AdS boundary. For us, it is easiest to state these conditions 
in coordinates which suggest a canonical choice for both these vectors. We write the AdS 
metric in spherical polar coordinates on where it takes the simple familiar form 

(6) QAdS = - (1 + dt^ + (l + dr^ + r’^dO^i , 

with the asymptotic boundary corresponding to the timelike hypersurface “r = oo” Note 
that 1/r is a boundary defining function and that 

eo = , ^ ^ dt , Cr = Vl -k r'^dr , eA = 

Vl + 


with [A = 1, 2) an orthonormal frame on the sphere of radius r defines an orthonormal 
frame for AdS. Finally, the vector dt singles out a preferred timelike direction. 

In the case of the wave equation, the optimally dissipative boundary condition can then 
be stated as 


(7) 


d{ru) 2 ^(™) 


-k r 


0 , 


as r 


oo. 


4.1 


dt dr 

Note that -k (l -k r^) dr is an outgoing null-vector. See also Section 

In the case of Maxwell’s equations, we define the electric and magnetic field Fj = 
F(eo,ei) and Hi = Ci) respectively. Here ^g^ds Hodge dual with re¬ 

spect to the AdS metric. The dissipative boundary condition then takes the form 


( 8 ) 


[Ea + €a^Hb) —> 0, as r 


oo, 


1 

2 


To simplify the algebra we choose A = —3 throughout the paper. 

The causal nature of this boundary is clearer in the Penrose picture discussed in Section 


1.5 


below. 
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which means that the Poynting vector points outwards, allowing energy to leave the space- 
time. 

Finally, in the case of the Bianchi equations, we define the electric and magnetic part of 
the Weyl tensor Eab = W (cq, ca, eo, es) and Hab = *gAds^ (eo, ga, eo, cb) respectively. 
The dissipative boundary conditions can then be expressed as 

(9) (eab ~ 2 ^abEc^ + j 0; as r —>■ oo . 


Interpreting the Bianchi equations as a linearisation of the full vacuum Einstein equations, 
we can understand these boundary conditions in terms of the metric perturbations as 
implying a relation between the Cotton-York tensor of the conformal metric on ^ and the 
“stress-energy tensor” of the boundar}]^ 

That the above boundary conditions are indeed correct, naturally dissipative, boundary 
conditions leading to a well-posed boundary initial value problem will be a result of the 
energy identity. While this is almost immediate in the case of the wave equation and 
Maxwell’s equations, we will spend a considerable amount of time on the “derivation” of 
in the Bianchi case, see Section [4.3[ 


We now turn to the results. 


1.2. The main theorems. 


Theorem 1.1. Let one of the following hold 

(W) is a scalar function and a smooth solution of ^ subject to dissipative boundary 
condition We associate with T the energy density 

£ [T] := \/l-F + dr {\/1 + r'^^^ -k 

(M) T is a Maxwell-two-form and a smooth solution of 0 subject to dissipative bound¬ 
ary conditions We associate with T the energy density 

e[^-] = Vl + r2(|E|2 + |Fr|2) 

where E and H denote the electric and magnetic part of T respectively. 

(B) T is a Weyl-field and a smooth solution of 0 subject to dissipative boundary con¬ 
ditions &■ We associate with T the energy density 

£[T] = {l + r'^Y (|E|2 + |F|2) 

where E and H denote the electric and magnetic part of the Weyl-field respectively. 


Then we have the following estimates 

(1) Uniform Boundedness: For any 0 < T < oo we have 


Vl + r 


-.r'^drdoj < 


£[T] 


Vi 


-.r’^drdui , 

I ■ 


where the implicit constant is independent ofT. 


^We use here the nomenclature of the putative AdS/CFT correspondence. 
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(2) Degenerate (near infinity) integrated decay without derivative loss: 


[ dt [ 

Jo J Yjf 


1 + 


r'^drdoj < 


L 


\/l + r' 


zr'^drduj. 


(3) Non-degenerate (near infinity) integrated decay with derivative loss: 


/o JUt 


Vl + r' 


-.r'^drduj < 




'So Vl + r 


Remark 1. Similar statements hold for higher order energies by commuting with dt and 
doing elliptic estimates. As this is standard we omit the details. 


Corollary 1.2 (Uniform decay). Under the assumptions of the previous theorem the fol¬ 
lowing uniform-in-time decay estimate holds for any integer n > 1 



r'^drduj < 


1 

(T+tf 



Vl + r2 


What is remarkable about the above theorem is that the derivative loss occurring in (3) 
allows one to achieve integrated decay of the energy without loss in the asymptotic weight 
r. While it is likely that more refined methods can reduce the loss of a full derivative in (3), 
we shall however establish that some loss is necessary and in fact reflects a fundamental 
property of the hyberbolic equations on AdS: the presence of trapping at infinity. 


Theorem 1.3. With the assumptions of Theorem [B the term s[dt^] on the right hand 
side of estimate (3) of Theorem \l.l\ is necessary: The estimate fails (for general solutions) 
if it is dropped. 


We will prove Theorem 1.3 only for the case of the wave equation (W), see (44) and 
Corollary 5.8[ The proof is based on the Gaussian beam approximation for the wave 
equation]^ In particular, we construct a solution of the conformal wave equation in AdS 
which contradicts estimate (3) of Theorem o without the e [St T]-term on the right hand 
side. Similar constructions can be given for the Maxwell and the Bianchi casej^ 


1.3. Overview of the proof of Theorem |1.1| and main difficulties. The proof of 
Theorem 1.1 is a straightforward application of the vectorfield method once certain diffi¬ 
culties have been overcome. Let us begin by discussing the proof in case of the wave- (W) 
and Maxwell’s equation (M) as it is conceptually easier. 

In the case of (W) and (M), statement (1) follows immediately from integration of the 


divergence identity V“ (Tab (dt) ) = 0 with T being the energy momentum tensor of the 


scalar or Maxwell field respectively. In addition, in view of the dissipative condition, this 
estimate gives control over certain derivatives of u (certain components of the Maxwell 
field, namely Ea and Ha) integrated along the boundary. 

The statement (2) of the main theorem is then obtained by constructing a vectorfield 
X (see ®) which is almost conformally Killing near infinity. The key observations are 


^See |17| for a discussion of the Gaussian beam approximation on general Lorentzian manifolds. 
^These will be more involved in view of the fact that the equations involve constraints. 
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that firstly, the right hand side of the associated divergence identity V“ = '’‘-vr • T 

controls all derivatives of u in the wave equation case (components of F in the Maxwell 
case). Secondly, when integrating this divergence identity, the terms appearing on the 
boundary at infinity come either with good signs (angular derivatives in the case of the 
wave equation, and Hr components in the Maxwell case) or are components already 
under control from the previous Tab (5t)^-estimate. The integrated decay estimate thus 
obtained comes with a natural degeneration in the r-weight, as manifest in the estimate 
(2) of Theorem o 

To remove this degeneration we first note that in view of the fact that the vectorfield dt 
is Killing, the estimates (1) and (2) also hold for the 5t-commuted equations. In the case 
of the wave equation, controlling dtdttp in on spacelike slices implies an estimate for all 
spatial derivatives of V’ through an elliptic estimate. The crucial point here is that weighted 
estimates are required. Similarly, one can write Maxwell’s equation as a three-dimensional 
div-curl-system with the time derivatives of E and H on the right hand side. Again the 
crucial point is that weighted elliptic estimates are needed to prove the desired results. 

Once all (spatial) derivatives are controlled in a weighted L?' sense on spacelike slices one 
can invoke Hardy inequalities to improve the weight in the lower order terms and remove 
the degeneration in the estimate (2). 

For the case of the spin 2 equations (B), the proof follows a similar structure. However, 
the divergence identity for the Bel-Robinson tensor (the analogue of the energy momentum 
tensor in cases (W) and (M)) alone will not generate the estimate (1). In fact, the term 
appearing on the boundary after integration does not have a sign unless one imposes 
an additional boundary condition! On the other hand, one can show (by proving energy 
estimates for a reduced system of equations, see Section 4.3.1) that the boundary conditions 
([^ already uniquely determine the solution]^ The resolution is that in the case of the 
Bianchi equations one needs to prove (1) and (2) at the same time: Contracting the Bel- 
Robinson tensor with a suitable combination of the vectorfields dt and X ensures that the 
boundary term on null-infinity does have a favorable sign and so does the spacetime-term 
in the interior. Once (1) and (2) are established, (3) follows from doing elliptic estimates 
for the reduced system of Bianchi equations similar to the Maxwell case (M). 


1.4. Remarks on Theorem ll.il The estimates of Theorem 11.11 remain true for a class of 
perturbations of AdS which preserve the general properties of the deformation tensor 
of the timelike vectorfield V^dt -\- X, cf. the proof of Proposition 5.6 This includes per¬ 


turbations which may be dynamical. This fact is of course key for the non-linear problem 


The estimates of Theorem 1.1 are stable towards perturbations of the optimally dis¬ 
sipative boundary conditions. In the cases of (W) and (M) one can in fact establish 
these estimates for any (however small) uniform dissipation at the boundary, cf. Section]^ 
Whether this is possible also in case of (B) is an open problem and (if true) will require 
a refinement of our techniques. In Section |6.2| we discuss the case of Dirichlet conditions 


®The introduction of this system in the context of the full non-linear problem goes back to [S]. 
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for (B) which fix the conformal class of the induced metric at infinity to linear order. The 
Dirichlet boundary conditions may be thought of as a limit of dissipative conditions in 
which the dissipation vanishes. We outline a proof of boundedness for solutions of the 
Bianchi equations in this setting. We also briefly discuss the relation of our work to the 
Teukolsky formalism in Section |6.3[ and argue that the proposed boundary conditions of 
m may not lead to a well posed dynamical problem. We state a set of boundary conditions 
for the Teukolsky formulation that does lead to a well posed dynamical problem. 


Finally, one may wonder whether and how the derivative loss in Theorem 1.1 manifests 
itself in the non-linear stability problem. As we will see in our forthcoming [18], the 
degeneration is sufficiently weak in the sense that the degenerate estimate (2) will be 
sufficient to deal with the non-linear error-terms. 


1.5. Main ideas for the proof of Theorem |1.3| , In order to better understand both 
the geometry of AdS and the derivative loss occurring in (3) of Theorem it is useful 
to invoke the conformal properties of AdS and the conformal wave equation Q. Setting 
r = tamp with pj G (O, |) one finds from ^ 

(10) gAdS = Ipdnl2) =: : 


cos- 


' p; 


COS- 




and hence that AdS is conformal to a part of the Einstein cylinder, namely M x C M x S^, 
where denotes a hemisphere of with V’ = f being its equatorial boundary. The wave 
equation Q is called conformal or conformally invariant because if u is a solution of Q, 
then V ;= flu is a solution of a wave equation with respect to the conformally transformed 
metric g = VP^gAdS- This suggests that understanding the dynamics of ([^ for gAdS is 
essentially equivalent to that of understanding solutions of 

( 11 ) Ug^v-v = t) 

on one hemisphere of the Einstein cylinder. This latter is a finite problem, which we will 
refer to as Problem 1 below. 

(1) Problem 1: The wave equation ( |ll[ ) on x (with the natural product metric 
of the Einstein cylinder) where is the (say northern) hemisphere of the 3-sphere 

with boundary at p) = |, where (say optimally) dissipative boundary conditions 
are imposed. We contrast this problem with 

(2) Problem 2; the wave equation on x (with the flat metric) where B^ is the 
unit ball with boundary where dissipative boundary conditions are imposed. 

For Problem 2 it is well-known that exponential decay of energy holds without loss of 
derivatives, see |20[ [211122j . (It is an entertaining exercise to prove this in a more robust 
fashion using the methods of this paper.) For Problem 1 , how ever, there will be a derivative 
loss present in any decay estimate, as seen in Theorem |l.lP| 

This phenomenon can be explained in the geometric optics approximation for the wave 
equation. Recall that in this picture, the optimally dissipative boundary condition says 


'It is straightforward to translate the estimates established for (j^ in Theorem 
picture currently discussed. 


1.1 


into the conformal 
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that the energy of a ray is fully absorbed if it hits the boundary orthogonally. For rays 
which graze the boundary, the fraction of the energy that is absorbed upon reflection 
depends on the glancing angle: the shallower the incident angle, the less energy is lost in 
the reflection. 

Now let us fix a (large) time interval [0, T] for both Problem 1 and Problem 2. To 
construct a solution which decays very slowly, we would like to identify rays which a) hit 
the boundary as little as possible and b) if they do hit the boundary, they should do this 
at a very shallow angle (grazing rays). 

What happens in Problem 2, is that the shallower one chooses the angle of the ray, the 
more reflections will take place in [0, T]. This is easily seen by looking at the projection of 
the null rays to the surface t = 0 (see figure). 




In sharp contrast, for Problem 1, the time until a null geodesic will meet the boundary 
again after reflection does not depend on the incident anglel This goes back to the fact 
that all geodesics emanating from point on the three sphere refocus at the antipodal point. 
As a consequence, in Problem 1 we can keep the number of reflections in [0, T] fixed while 
choosing the incident angle as small as we like. This observation is at the heart of the 


Gaussian beam approximation invoked to prove Theorem 1.3 


1.6. Structure of the paper. We conclude this introduction providing the structure of 
the paper. In Section we define the coordinate systems, frames and basic vectorfields 
which we are going to employ. Section introduces the field equations of spin 0, 1 and 2 
fields, together with their energy momentum tensors. The well-posedness under dissipative 
boundary conditions for each of these equations is discussed in Section with particular 
emphasis on the importance of the reduced system in the spin 2 case. Section is at the 
heart of the paper proving the global results of Theorem [Hi first for the spin 0 (Section 
5.1), then the spin 1 (Section |5.2[) and finally the spin 2 case (Section |5.3[). Corollary 1.2 


is proven in Section [5.4| and Theorem 1.3 in Section [5.5[ We conclude the paper outlining 
generalizations of our result. Some elementary computations have been relegated to the 
appendix. 
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2. PRELIMINARIES 

2.1. Coordinates, Frames and Volume forms. We will consider the so-called “global” 
anti-de Sitter space-time, defined on with metric (in polar coordinates) 

dr ^ 

(12) gAds = “ (l + ^ ^2 

with d ^'^2 the standard round metric on the unit sphere. As long as there is no risk 
of confusion, we will also denote the metric by g. It will be convenient to introduce an 
orthonormal basis {e“}: 

(13) e° = VT+T^dt, e” = , , 

yjl + r'^ 

for A = 1,2, where are an orthonormal basis for the round spher^of radius r. Through¬ 
out this paper, we will use capital Latin letters for indices on the sphere while small Latin 
letters are reserved as spacetime indices. The dual basis of vector fields is denoted {cq}: 

1 ^ ^ 

(14) = 

We introduce the surfaces = {t = T} and = {r = R,Ti < t < T 2 }, which 

have respective unit normals: 

n := eo, rn := e^. 

The surface measures on these surfaces are 

^2 _ 

d5sj, = , drdu!, dS-^n = v 1 -|- r'^dtduj 

with do; the volume form of the round unit sphere. We denote by ■= {Ti < t < T 2 } 

the spacetime slab between T,Ti and Finally, the spacetime volume form is 

dg = r^dtdrdoj. 

2.2. Vectorfields. The global anti-de Sitter spacetime enjoys the property of being static. 
The Killing held 

(15) T := — = \/r+r^eo 

is everywhere timelike, and orthogonal to the surfaces of constant t. 

Besides the Killing held T, we will exploit the properties of the vectorheld 

(16) X := r\/l + r"^— = rcr 


O o 

Obviously there’s no global orthonormal basis for S . We can either take multiple patches, or else 
understand the capital Latin indices as abstract indices. 
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While the vectorfield X is not Killing, it is almost conformally Killing near infinity. More 
importantly, it generates terms with a definite sign. To see this, note that 

-I _ 

(17) ^TT:=-Cxg= , - + ff\/l + r2. 

^ Vl + T 

In particular, contracting ^vr with a symmetric traceless tensor T will only see the first term 
and this contraction will in fact have a sign if T satisfies the dominant energy condition. 

2.3. Differential operators. Throughout this paper, we will use V to denote the Levi- 
Civita connection of g. Define also as the standard Laplace-Beltrami operator: 


□gU := V“VaM. 


2.4. The divergence theorem. We denote by Div the spacetime divergence associated 
to the metric g for a vector field. In coordinates, it takes the following form: 

If K = KfOeo + then 

(18) DivK = -^j^^ + lT(PvTT75/U)+'!:V4. 

It will be useful to record the following form of the divergence theorem: 

Lemma 2.1 (Divergence theorem). Suppose that K is a suitably regular vector field defined 
on M^. Then we have 


0 = 


'St, 


Kan°-dS^^^ - 


KaU'^dST. 


' Stj 


Tl 


— lim / 

r^oo 


KamSdSf^^ + 


Div Kdg. 


dTi.Tal 


3. The field equations 

3.1. Spin 0: The wave equation. The conformal wave equation on AdS is 


(19) 


□ 


pAdS 


u + 2u = 0. 


A key object to study the dynamics of (19) is the twisted, or renormalised energy-momentum 
tensor. See [231 for further details on the following construction. We first define the 
twisted covariant derivativ^ by: 

With this derivative we construct the twisted energy momentum tensor 


( 20 ) 


1 


TafeM = VaUVbU - -gab VcUVu + 


1 -I- 


'^When acting on a scalar, we will sometimes write V^/ =: 5^/. 
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which satisfies 

(21) = {n\gU + 2u) VbU - T'^cVbl. 

If K is any vector field then we can define the current 

(22) ^r[n] =T^[n]i^^ 

which is a compatible current in the sense of Christodoulou [25]. Moreover, if K is Killing 
and K{r) = 0, then = 0 and we can see that is a conserved current when u 

solves the conformal wave equation (19). 

3.2. Spin 1: Maxwell’s equations. The Maxwell equations in vacuum are given by the 
following first order differential equations for a 2-form F on 

dF = 0, d {-*:gF) = 0. 

We decompose the Maxwell 2—form as: 

(23) F = Kj. e° A + Ka e° A + Hr A + Ha A e^, 

where cab is the volume form on the spherj^ The dual Maxwell 2—form is: 

*gF = Hre° A e^ + Ha e° A - Er A - Ea A e^. 

We will use the notation Ei with i = f, 1, 2, and similarly for H. Smoothness implies that 
Ei, Hi are bounded (but not necessarily continuous) at the origin. With respect to this 
decomposition, Maxwell’s vacuum equations (3.2) split into six evolution equations: 

(24) -j=L=dtEr = -re^^l/AHB, 


(25) 

(26) 
(27) 


\/rT 

r 

r 

■y/l +1 


-AEa = eA 


B 


dr (t^I^Hb) -T^bH, 


-AHr = 


re 


AB 


A^B-, 


\/l + T- 


AH A = -eA 


B 


dr {r \/1 -I- t'^Eb] — r'f B^i 


and two constraints: 
(28) 

(29) 


0 = 


0 = 


Vl + r'^ 

r 

Vl + r"^ 


dr {r^Er) + rf^EA, 

dr {r^Hr) + rf^HA. 


Here is the covariant derivative on the unit sphere, which commutes with The 

evolution equations (24 - 27) form a symmetric hyperbolic system. If the evolution equations 


hold, and assuming sufficient differentiability, it is straightforward to verify that 

d^_dJ^_ 

At ~ 


10 


ei 2 = 1 in this basis. 
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where 

•= ^^^''^ dr{r^Er)+rf^EA, 

Jif := ^ll±^dr {r^Hr) + rf^HA. 

Thus if the constraint equations hold on the initial data surface, then they hold for all 
times. 


3.2.1. The energy momentum tensor. The analogue of the energy momentum tensor (20) 
for the wave equation is the symmetric tensor 


(30) 


nPU = FacEb^ - IgabFcdF^^, 


which satisfies 

V“T[F],fe = + {*F)b‘^V\*F)dc 

so that if F satisfies Maxwell’s equations, T[F] is divergence free and traceless. We define 
in the obvious fashion the current 


(31) ^r[F] = T'^b[F]K^ 

for any vector field K. 


3.3. Spin 2: The Bianchi equations. The equations for a spin 2 field, also called the 
Bianchi equations, can be expressed as first order differential equations for a Weyl tensor, 
which is an arbitrary 4—tensor which satisfies the same symmetry properties as the Weyl 
curvature tensor. More precisely: 


Definition 1. We say a 4:—tensor W is a Weyl tensor if it satisfies: 
i) WaJicd — bffeacd — ^^abdc- 

ii) abed T Wacdfe + W^abde — 0. 

Hi) W^abed — fTcciaft- 
iv) W°-bad = 0. 

The dual of a Weyl tensor is defined by (ef. [I]J 

^^abed ^ fed ■ 

The Bianchi equations are 
(32) V[aWbc]de = 0, 

which are equivalent to either of the equations 


(33) 

(34) 


y^Wabed = 0, 
y"*Wabcd = 0 , 


and for us studying (33) will be particularly convenient. 

As in the case of the Maxwell equations, it is convenient to decompose the Weyl tensor 
based on the 3 + 1 splitting of space and time. In the Maxwell case, the field strength 
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tensor F decomposes as a pair of vectors tangent to In the spin 2 case, the Weyl tensor 
W decomposes as pair of symmetric tensors tangent to We dehne: 


(35) Eab := WoaOb = hhcadfe(eo)'^(eo)'^, 

(36) Hab := *WoaOb = ^WcadbieoTieof. 


The symmetries indeed ensure that E, H are symmetric and tangent to S*. Moreover, both 
E and H are necessarily trace-free. In fact, one can reconstruct the whole tensor W from 
the symmetric trace-free fields E, El, see [U §7.2]. 

We will further decompose E, H along the orthonormal frame defined in ([T^. To write 
the equations of motion in terms of E, H we consider the equations 0 = V^'Wawo, 0 = 
V“Wa(Ar)0 and 0 = V“IT„(Aij)o, from which we respectively find the evolution equations 
for E: 


(Evol Ej^) 
(Evol Eat) 

(Evol Eab) 


r dErf 

Vl +1-2 dt 

r BEat 

Vl + dt 

r BEab 
Vr+72 dt 


= -re^^fAHBf, 

1 B 

dr {r{l + r‘^)HBf) 

^/l + r'^ 

- 

dr (r(l -h r‘^)HB)c) 

e(A 

Vl + 


B^rr 




From the equations 0 = V“ITaOAo and 0 = V“ITaOrO respectively we find the constraint 
equations: 

(Con Er) ^ if Erf) + t^^Ebt =: = 0 

(Con Ea) ^ + rf^EAB =: A = 0 

By considering the equivalent equations for *W, we find that the evolution equations for 
H can be obtained from these by the substitution {E,H) —)• {H, —E): 


(Evol Hff) 

r dHff 

= re^^fAEBr, 

Vl + dt 

(Evol HAr) 

r dHAr 

1 B 

dr (r(l -h r‘^)EBf) 

Vl + 7-2 dt 

Vl + 

(Evol Hab) 

r BHab 

- 

dr (r(l r‘^)EB)c) 

Vl -h r2 dt 

Ha 

Vl + 


T ^ bEtt 

- r^\c\EBy 


+ "^e^^fBEcA, 
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and the constraint equations for H are: 

(Con Hr) ^ {r^Hrr) + rl/^HBr =: ^ = 0, 

(Con Hb) {r^HAr) + tU^Hab =: = 0. 

3.3.1. The Bel-Robinson tensor. The spin 2 analogue of the energy momentum tensor for 
the Maxwell field is the Bel-Robinson tensor [H §7.1.1]. This is defined to be 

(37) Qabcd := WaecfW, ^ + ^Waecf^W, ^ 

It is symmetric, trace-free on all pairs of indices and if W satisfies the Bianchi equations 
then 

(38) V^Qabcd = 0. 

We shall require the following quantitative version of the dominant energy condition for 
the Bel-Robinson tensor 


Lemma 3.1. Suppose that ti,t 2 G TpM. are future-directed timelike unit vectors, and that 
—g{ti,t 2 ) < B for some B >1. Then there exists a constant C > 0 such that 

■^^Q{ti,ti,ti,ti) < Q{ti,tj,tk,ti) < CB'^Q{ti,ti,ti,ti) 

holds for any i,j,k,l G {1,2}. Moreover, Q{ti,ti,ti,ti) > 0, with equality at a point p if 
and only if W vanishes at p. 


Proof. We follow the proof of m Prop. 4.2]. We can pick a pair of null vectors e§, 64 with 
g{e'^,e'^) = such that ti = e '3 -|- 64 and t 2 = be'^ b~^e'^ for some 6 > 1. The condition 
— 5 (^ 1 , ^ 2 ) < B implies 

i<.<R 

As a consequence, we can write 


Q{ti,tj,tk,ti) = gr“/fcfQ(e'„,e{,e(,,erf) 


a,6,c,d=3 


where 

^ ^ „abcd ^ / rj4 

JfA ^ ^ijkl <CB 

for some combinatorial constants c,c' which are independent of B. By |16[ Prop. 4.2], 
all of the quantities Q(e(j, e}, e),, e'^) are non-negative, and so we have established the hrst 
part. For the second part, we note that, again following m Prop. 4.2], the quantity 
Q{ti,ti,ti,ti) controls all components of W. □ 


4. Well posedness 


4.1. Dissipative boundary conditions. 







ASYMPTOTIC PROPERTIES OF LINEAR FIELD EQUATIONS IN ANTI-DE SITTER SPACE 


15 


4.1.1. Wave equation. For the wave equation, an analysis of solutions near to the conformal 
boundary leads us to expect that u has the following behaviour: 


u+{t,x^) , u-{t,x^) , I 
u =-1-^- \- O I ^ 


as r —>■ oo. 


The boundary condition we shall impose for (19) on anti-de Sitter space can be simply 
written in the form 


(39) 

which is equivalent to 


d{ru) o djru) 
dt dr 

du. 


as r —> oo, 


dt 


— n_ = 0 , on . 


4.1.2. Maxwell’s equations. In the Maxwell case, by (optimally) dissipative boundary con¬ 
ditions we understand the condition 

(40) (^Ea + eA^Hs) —>-0, as r —oo . 

The above boundary conditions ensure that the asymptotic Poynting vector is outward 
directed, i.e. energy is leaving the spacetime. It effectively means that behaves like 
an imperfect conductor with a surface resistance, for which the electric currents induced 
by the electromagnetic field dissipate energy as heat. These boundary conditions are an 
example of Leontovic boundary conditions [26( §87]. 


4.1.3. Bianchi equations. In the case of the Bianchi equation, by (optimally) dissipative 
boundary conditions we understand the condition 

(41) ^Eab ~ 2 ^abEc^ + 0; as r —>■ oo . 

Notice that there are only two boundary conditions imposed (since the relevant objects 
appearing in the boundary term are the traee-free parts of Eab, Hab)- We could choose 
more general boundary conditions, for the time being we shall just consider those above, 
which ought in some sense to represent “optimal dissipation” at the boundary. 


4.2. The well-posedness statements. We now state a general well-posedness statement 
for each of our three models with dissipative boundary conditions. As is well-known, the key 
to prove these theorems is the existence of a suitable energy estimate under the boundary 
conditions imposed. In the Wave- and Maxwell case such an estimate is immediate. In 
the Bianchi case, however, there is a subtlety (discussed and resolved already in | 8 ]). We 
will dedicate Section 4.3 to derive a local energy estimate showing that the condition (41) 
indeed leads to a well-posed problem. 
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4.2.1. Spin 0: The wave equation. The following result can be established either directly, 
or by making use of the conformal invariance of (19); 


Theorem 4.1 (Well posedness for the conformal wave equation). Fix T > 0. Given smooth 
functions uo,ui : Sq —)> M satisfying suitable asymptotic conditions, there exists a unique 
smooth u, such that: 


i) u solves the conformal wave equation (19) in 
a) We have the estimate: 


sup ^ 

*[0T] k,l,m<K 


(dif {rf)"^ (ru) < Ct,uo,ui,k 


for k,l,m = 0,1,... and a constant Ct,uo,ui,k depending on K, T and the initial data. 
This in particular implies a particular asymptotic behaviour for the fields. 

Hi) The initial conditions hold: 


u\ 


t=0 ~ ^0) 




iv) Dissipative boundary conditions (39) hold. 


4.2.2. Spin 1: Maxwell’s equations. Working either directly with (24- 29), or else by making 
use of the conformal invariance of Maxwell’s equations, it can be shown that: 


Theorem 4.2 (Well posedness for Maxwell’s equations). Fix T > 0. Given smooth vector 
fields E^,H^ satisfying the constraint equations, together with suitable asymptotic condi¬ 
tions, there exists a unique set of smooth vector fields: Ei{t),Hi{t), such that: 


i) Ei{t), Hi{t) solve (24~ 29) in 5'[o^t] corresponding Maxwell tensor E is a smooth 

2—form on 

a) We have the estimate: 


sup 

S[0,T] 


kl.m ^ j. 




{rUyir^Ef) + {r^drf {dfj' irUY' [r^Ef) 




for any K > 0, where Crp^E0,H0,K depends on K,T and the initial data, 
particular implies a particular asymptotic behaviour for the fields. 

Hi) The initial conditions hold: 


Em=El Hm = H^ 


iv) Dissipative boundary conditions (40) hold. 


This in 


The asymptotic conditions on the initial data are corner conditions that come from 
ensuring that the initial data are compatible with the boundary conditions. It is certainly 
possible to construct initial data satisfying the constraints and the corner conditions to 
any order. We could work at finite regularity, and our results will in fact be valid with 
much weaker assumptions on the solutions, but for convenience it is simpler to assume the 
solutions are smooth. 
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4.2.3. Spin2 : Bianchi equations. In the case of the Bianchi equation we can prove: 


Theorem 4.3 (Well posedness for the Bianchi system). Fix T > 0. Given smooth trace¬ 
less symmetric 2-tensors on Sq satisfying the constraint equations, together with 

suitable asymptotic conditions as r —)■ oo, there exists a unique set of traceless symmetric 
2-tensors: Eab{t),Hab{t) such that: 

i) Eab{t), Hab{t) are traceless and the corresponding Weyl tensor W is a smooth 4—tensor 
satisfying the Bianchi equations on 5'jo^'r]. 
a) We have the asymptotic behaviour: 


sup ^ {r‘^dr)^ {dt^ {rf)'^ {r^Eab) + ( 

'S'[o,T] 


{dtY {rUr' {r^H^b) 


< CrpE0^H0,K 


k'V 


,<K 


as r ^ oo, for any K >0. 

Hi) The initial conditions hold: 

Eab{0) = E^^, 


Hab{0) = H\ 


0 

ab 


iv) Dissipative boundary conditions (fl) hold. 

We will spend the remain der of this section to derive the key-energy estimate that is 
behind the proof of Theorem 4.3^^ 


4.3. The modified system of Bianchi equations. In order to establish a well posedness 
theorem for the initial-boundary value problem associated to the spin 2 equations, the 
natural thing to do is to consider the evolution equations (Evol) as a symmetric hyperbolicp^ 
system. Having established existence and uniqueness for this system (Step 1), one can then 
attempt to show that the constraints (Con) are propagated from the initial data (Step 2). 

If one attempts this strategy with the equations in the form given in Section |3.3t one 
finds that Step 1 causes no problems: one can easily derive an energy estimate for the 
system (Evol). In fact, it is a matter of simple calculation to check that taking: 


r(l -|- r ^)2 X (Evol Ej^) + 2EAf x (Evol Eat) + Eab x (Evol Eab ) + E H] 
and integrating over with measure drduj we arrive alp^ 

d 1 r 

dt 2 Js, 


-|- 2 E r^(l -|- r‘^)drduj 


(42) 


= lim / (\e^^EArHBr + (-^^EacHb^^ r^div . 


^ ^Theorem of course establishes a stronger (global) estimate. The key point here is to explain why 
the naive approach using only the Bel-Robinson tensor and the Killing fields fails and also to der ive the 
(reduced) system of equations which will be needed in the second part of the proof of Theorem |l.l[ 

^^Or at least Friedrichs symmetrizable, but the distinction is not important for our purposes. The key 
issue is the existence of a good energy estimate. Providing this exists, putting the system into symmetric 
hyperbolic form is straightforward. 

^^Here and elsewhere, in expressions like \EAr\'^ -f \Eab\^, contraction over the indices A,B etc. is 
implied, i.e. -b {Eab^ = EArE^r + EabE^^ 
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Once one has an energy estimate of this kind, establishing the existence of solutions to 
the evolution equations under the assumption that the boundary term has a good sign is 
essentially straightforward. Looking at the boundary term we obtain, this suggests that 
boundary conditions should be imposed on both the Ar and the AB components of E (or 
alternatively H) at infinity. We shouldn’t be so hasty, however. Before declaring victory, 
we must return to look at the constraints (Step 2). 

Firstly, it is simple to verify that if (Evol) hold with sufficient regularity then 

_ dEa'^ _ dHa^ 

dt dt ’ 


so that the trace constraints on E and H are propagated by the evolution equations. Next, 
we turn to the differential constraints. We find that if (Evol) and the trace constraints 
hold, then the functions S‘a,J^a defined in Section [3(3| satisfy the system of equations; 


(Evol Sp) 

r dSp 

_ ^ 




\/l + r2 dt 

2 




(Evol Sa) 

r SSa 


dr (^r^Vl + 

- - 

JCB 

Vl + r’2 dt 

2 


\/l + 

\/l + 

(Evol .Jifp) 

r dJ^ 

_ L^AB 

^aS’b , 



.y/l _|_ 7-2 Dt 

2 



(Evol JSa) 

r dJ^A 


dr + t'^Sb^ 

B'^r 

Sb 

Vl + r2 dt 

2 


\/l + 



Now, things appear to be working in our favour. This system is symmetric hyperbolic and 
we can check that by taking: 

r^(l + r^) [Sr X (Evol Sp) + Sa x (Evol Sa) + x (Evol + ASa x (Evol JSa)] 
and integrating over with the measure drdoj we can derive: 
d 1 


dt 2 


L 


ISrl"^ + ISa^ + + \JifAf r^yjl + r‘^drduj 


= lim 
r^oo 


-e 


r^duj— / [e^^S a^b\ r^\/l + r'^drdu:. 

's^nEt / Jilt 

Here we see the problem. We have no reason a priori to expect that the boundary term 
on Tjr n Tit vanishes. If it did, we could infer by Gronwall’s lemma that the constraints 


are propagated. We conclude that the form of the propagation equations of Section 3.3 


does not in general propagate the constraints at the boundary if boundary conditions are 
imposed on both Eab and Eap (or Hab and Hap)- 


4.3.1. The modified equations. In order to resolve this issue, we have to modify the propa¬ 
gation equations before attempting to solve them as a symmetric hyperbolic system. In the 
previous calculation, the problematic boundary terms arise due to the radial derivatives 
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appearing on the right hand side of (Evol Eat), (Evol Hat)- We can remove these radial 
derivatives at the expense of introducing angular derivatives by using the constraint equa¬ 
tions (Con Ejr). It is also convenient to eliminate Eyf and Hrf from our equations using 
the trace constraints. Doing this, we arrive at the modified set of propagation equations; 

eA^HBf 


(Evol’ Eat) 
(Evol’ Eab) 
(Evol’ Hav) 
(Evol’ Hab) 


r 

OEat 

Vl + 

! dt 

r 

OEab 

Vl + 

dt 

r 

dHAr 

Vl + 

dt 

r 

dHAB 

Vl + 

dt 


= —re 


BC 


bHcA - 


VTTr^’ 


= e 


c 


{A 


1 


d 


Vl + dr 

= re^^fs^CA + 

1 d 


(r(l -b r‘^)HB)c) - r^\C\HB)f 
^A^Ebt 


Vl + 


= —e 


c 


{A 


(r(l -h r^)EB)c) - rf\c\EB)f 


Vl + dr 

This again forms a symmetric hyperbolic system. Taking 

r(l -|- r ^)2 [Eav X (Evol Eat) T Eab ^ (Evol Eab ) T E -fA id] 
and integrating over with measure drdu we arrive at 


^1 r 
dt 2 J-£ 


E\Eab\ E\HAr\ E\Hab\ r"^^! + r‘^)drdio 


(43) 


= lim 

r^oo 


^EAcHB^r^duj + 2 


s^nSt 


^EAfHBr l(1 + r‘^)drduj , 


which is certainly sufficient to establish a well posedness result for the equations (Evol’), 
provided we choose boundary conditions such that the term on has a sign. Notice 
that this will involve imposing conditions only on the Eab (or Hab) components, so this 
formulation of the propagation equations is clearly different to the previous onep | 

Notice also that, unlike the estimate (42) for the (Evol) equations, we now have a bulk 
term in the energy estimate (43). For well-posedness this is no significant obstacle, but it 
will make establishing global decay estimates more difficult. In particular, it is no longer 
immediate that solutions with Dirichlet boundary conditions remain uniformly bounded 
globally. 

Now let us consider the propagation of the constraints. We of course have to interpret 
the E^ term in V as {—Ea^) and similarly for The evolution equations for the 


14t 


■‘For instance, specifying Eab = 0 on J*’’ will clearly lead to a unique solution of the modified system. 
On the other hand, the same boundary condition will not completely fix the boundary term occurring in 
(421 in the unmodified formulation, illustrating the severe drawback of the latter. 
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constraints take a simpler form: 
(EvoP 4 .) 

(EvoP 44 ) 

(EvoP Mf) 

(EvoP 


94 


\/l + r2 dt 
r dfoA 

VT +72 dt 

r dJ^ 


B 










= re 


,A_B 


41+72 

r 9+^ eA^ ’’ 




i?4 + 




VT +72 


47+72 Qt 2 

Now, once again this is a symmetric hyperbolic system. Taking 
r^{l + r^) [4 X dEvol’ 4D + 244 x dEvoP Sa\ + ^ x dEvol’ +4D + 2+4i x dEvoP 


and integrating over with the measure drdoj we can derive: 
d 1 


dt 2 
= -2 


|4|V 2 |44p + 1941^2 rVl +4drda; 

S^a^^b] r^\/l + r'^drduj 


Immediately, with Gronwall’s Lemma, we deduce that the constraints, if initially satisfied, 
will be satisfied for all time. 

To our knowledge, identifying the above modified system as the correct formulation to 


prove well-posedness goes back to Friedrich’s work [8]. In particular, Theorem 4.3 above 
could be inferred from this paper. 

5. Proof of the main theorems 


5.1. Proof of Theorem |1.1| for Spin 0. The Killing field T immediately gives us a 
boundedness estimate: 


Proposition 5.1 (Boundedness of energy). Let u be a solution of (19) subject to dissipative 
boundary conditions (39) as in Theoreni\4.1\ Define the energy to be: 


(44) 


Etlu] := - 


{dtu)^ + 


'St 


1 + ^2 


+ (1 + r^) (^drU^ + \'fu^ j r'^drdu. 


Then we have for any Ti < T 2 : 


EtM + ^ [r\dtu)^ + 4(44^ ) dujdt = EtM- 

'J Soo 


Proof. We have that 


Div (^J) = 0. 
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\'T' 'T' 1 

Integrating this over Pi^k up terms from and ^ . A straightforward 

calculation shows 





Et[u] 


We also find 


= f r^{l + r^) (dtu) {dru)dujdt 

u S-p J 

\^{dtuf+ r^{l+r^f{druf 
|5t(™) + r(l + r^)(9r.u)| 


if 

^ JS2 


dbjdt. 


As r —)> oo, the term in braces vanishes by the boundary condition and we find 


lim [ '^Iam°‘dSy = —- 

r^oo ylTl,T2] “ W 2 


/-[T1.T2] (r'^idtuf + r^{druf) dudt. 
Soo 


Applying Lemma 2.1, we are done. 


□ 


We next show an integrated decay estimate with a loss in the weight at infinity: 

Proposition 5.2 (Integrate decay estimate with loss). Let T 2 > Ti. Suppose u is a solution 
of (19) subject to as in Theorem j.l. Then the estimate 


{dtuf + 


'S[Ti,T2] 


1 _j_ J.2 


VT+~r 


zdrdudt 


+ (1 + r^) (^dru'^ + \ ^u 

+ f (r‘^{dtu)‘^ + \r‘^'^uf + r^{dru)‘^')dujdt<CETAu] 
/f.[T'l,T'2l V ' ' / 

O Zjoo 


holds for some constant C > 0, independent of Ti and T 2 . 

Proof. We integrate a current constructed from the (renormalized) energy momentum ten¬ 


sor (20) and a radial vector field. The current is 


ja 

^ X,Wi ,W 2 


= uV^u + W2U^X'^ 

V1 -|- 


where X is the radial vector field defined in (16). The proof of the theorem is a straight 


forward corollary of the following two Lemmas. 
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Lemma 5.1. We have 

Div Jx,wi,W 2 = (^J • eo) - - \/l + + 

+ (Va 


Wl 


\/rT^ 


T'" 


Wl 


+ {l + r^)Div 


+ 2w2Xa uV U 


W2 
1 + 


X - 


Wl 


(1 + 7-2)2 


u^. 


In particular, if we take wi = 1, W 2 = + r'^) ^ then we have 




Div Jx,wi,w 2 — ( J ■ ^o) “1“ ^ 

2(1 + r2) 2 

Proof. The vector field X has the deformation tensor; 


+0+ ^_ 

^TT = - + 5(\/l + r2, 


y/l+r 


so that 


Div • eo) - (x'’Va - Vl + r^) T'^c 

We also require the observation that 

/ ~ \ 

VaU 


\/l + r2V“ 


\/l + r2 


— (CHg-'U + 2u) + 


Thus when u solves the conformal wave equation, we have: 


Div 


Wl 


\/l + r 


= - 




\/l + r 




Wl 


\/l + r2 


1 + J.2 


uV“u - 


U^Wl 

(1 + r2)i 


Finally, we have 


Div {w 2 U^X) = 2w2uX^X aU + (1 + r^)Div 


W2 

1 _|_ j,2 


X u\ 


Combining these identities we have the first part of the result. We can arrange that the term 
proportional to T'^c vanishes by taking wi = \. The term proportional to uXaU vanishes if 
W 2 = 5 (l + r^)“^, and the final part of the result follows from a brief calculation. □ 

Lemma 5.2. With wi, W 2 chosen as in the second part of the previous Lemma, we have: 


I ’^X,Wi,W2^0,^^^T2 


< CEtM < CEtM 


and 


lim / Jx 1 IU mn'niadSf. + - / 

r-;.oo 7g[7’i.r2l ^’^i’-!«2 “ S,. 2 7g[7’i,7’2] 


r'^'fu dcudi < CFIti M- 
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Proof. We calculate 


dtu 


'st (1 + r2)2 


r(l + r^)drU + u ) r^drdoj 


which, after applying Cauchy-Schwarz, can certainly be controlled by the energy Et[u] 


which in turn is controlled by Eo[u\ using Theorem 5.1 
For the other surface terms, we calculate 


1 


J X,W\,W2''^O.^^T,r ~ Q 


/V'L 


2 /v[ri.r2] 


r'^idtu) , 2n , 2 


1 _j_ ^2 


+ r"'(l + r^){dru) + 


+ 2ru[dru) — Ir^r 


r\/l + r^dtdio 


so that 


lim / 

r^oa J^[ri,T2] 


'^X,Wl,W2^0,^‘^Er ~ 


1 


r‘^{dtu)‘^ — + r^{dru )‘^) dujdt 


2 E\Ti,T2] 

^ ■^OO 

and the result follows since we already control the time and radial derivatives of u on the 
boundary by Theorem 5.1 □ 


This concludes the proof of Proposition 5.2, and establishes the claimed degenerate 
integrated decay without derivative loss result. □ 

We next improve the radial weight of the spacetime term in the integrated decay esti¬ 
mate, at the expense of losing a derivative 


Proposition 5.3 (Higher order estimates). Let u be a solution of (19) subjeet to (39) as 
in Theorem\4.1\ Then the estimate 




(1 + r^)2 


[Ti,T2] 


dr ( 


u 


12 'f 

+ 


l + r^ 


dr [r^u] + dr [dtu] 


drdudt 


+ 


L 


\^‘^u\‘^drdijjdt < C{Et^ [n] + Et^ [ut]) 


^S[Ti,T2] ^ 

holds for some constant C > 0, independent of Ti and T 2 . 


Proof. Let us define, for a solution of the conformal wave equation: 

Lu := {utt + u) = dr (r'^dr (^u\/l + . 

By commuting with T and applying Proposition |5.2[ we have 

(45) f I [Luf \dr]<C {Et^ [m] + Et, [ut]) . 

J ■j’2] ^ V -L “r ^ J 

shall take the frugal approach of commuting with the timelike Killing field. If one is happy to 
exploit the angular symmetries of AdS, our approach can be simplified. 
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We can expand the integrand to give 


[Lu] 


VT+~r2 


VTTr^ 


Vl + 


- 2 

dr ^r^dr (u \/1 + 


(46) 


+ 2 


dr (r^dr (u^l + r^ 




We clearly have two terms with a good sign and a cross term. To deal with the cross term, 
we integrate by parts twice, so that we obtain a term (with a good sign) that looks like 
and some lower order terms. More explicitly, we have 

Lemma 5.3. Let K he the vector field given by 
K = dr (r^dr (^u ^/1 + ca — {r')l7j^u)dr 


U Cr — 


2(1 + r2) 


\r'^u\ Cr- 


Then we have 


dr [^r^dr (u\/l + = \/l + dr ir^u) 


+ 


3 Ir^ 


u\ 


2(l + r2)5 


+ Div K. 


Proof. See Appendix |7.1 


□ 


To prove Proposition |5.3[ we simply insert (46) into (45) and handle the cross-term nsing 
Lemma |5.3[ The boundary terms coming from Div K are 


Div Kdrj = — 


’[T1T2I 


IL 


[ri.Tj] 


duidt 


which we control by the estimate of Proposition 5.2 

This in particular controls the term . Since has constant positive Gauss curva¬ 

ture, we have the following elliptic estimates (see for instance Corollary 2.2.2.1 in HI): 

JS2 JS2 

from which we obtain the desired bounds for _ 

Finally, the {drdtu)'^ term appearing in Proposition 5.3 is directly controlled by the 
T-commuted version of the estimates in Proposition |5.2[ 

□ 


We finally improve the weigh t in the spacetime term of Proposition |5.2| making use of 
the fact that by Proposition ^ we now control radial derivatives of dtu, and drU which 
lead to improved zeroth order terms through a Hardy inequality. This is a standard result, 
but for convenience we include here a proof. 



































ASYMPTOTIC PROPERTIES OF LINEAR FIELD EQUATIONS IN ANTI-DE SITTER SPACE 


25 


Lemma 5.4 (Hardy’s inequality). Fix a ^ 0. Suppose that f : [l,oo) —)> M is smooth, 
/(I) = 0 and |/(t)|^ —>-0 as r —)• oo. Then we have the estimate 

/ OO A poo 

provided the right hand side is finite. 

Proof. We write 



Here we have used /(I) = 0 and the fact that linir-^-oo t “ |/(r)p = 0 to discard the 
boundary terms. Now applying Cauchy-Schwarz, we deduce 

/ OO /a poo poo \ ^ 

|/|2r-i-“dr< \drff r^-^drj , 

whence the result follows. □ 


From Lemma 15.41 we establish: 

Theorem 5.1 (Full integrated decay). Letu he a smooth funetion sueh that \ru\ is bounded. 
Then the estimate 


“'>S’[Ti,T2] 

< c 

+ 


dr (j’^dru'j 


^ ' r'^drdivdt 


il + r^)-2 


’[ri.Tal 


1 2 r 

+ ^ 
1 _|_ j,2 


(dtu)'^ + 
1 _|_ ^2 


<C{ETAu]+ETAdtu]). 


dr \r'fu 

+ (l + r^) (^rV^ + 


+ 


dr [dtu] 


drdojdt 


\/l + r2 


drdivdt 


holds for some constant C > 0 independent ofTi and T 2 . 


Proof. By introducing a cut-off we can quickly reduce to showing that the estimate holds 
for u supported either on r < 2 or on r > 1. For u supported on r < 2, the estimate follows 
immediately, since the first order terms on the right hand side are comparable to those on 
the left hand side on any finite region. For u supported on r > 1, we first apply Lemma 
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5.4 with / = r^Vl + r'^drU and a = 1 to deduce 


'5[TlT2l 


5, 


u 


r^(l + r^)drdujdt < 4 


< 4 


' 5 [Ti,T 2 ] 


'S[ 


[ 3 ’iT 2 ] 


dr{r ‘^\/1 + r^dru) 
(1 + r^) (9r 


drdudt 


1 2 


< 4 


f (l-hr^)2 

to 

'^ITi,T2] 

L \ / J 


drdudt 

drdcodt. 


Similarly, applying Lemma 5^ to / = \/l + r‘^dtu with a = 1 we deduce 

[ \dtu\^ 

J S'r'T. 'T.,! 


,2 


’[ri,T2] 


1 _|_ j,2 


drdiodt < 


\dtu\‘^ r + r^)drduj 


'S[Ti,T2] 


< 4 


< 


L 


[T 1 .T 21 


'S[Ti,T2] 


dr[dtu\ (1 + r'^)drdujdt 
2 , - 

dr[dtu] r^yl + r'^drdujdt, 


where we’ve used that u is supported on r > 1. A similar calculation gives the estimate 
for yu. □ 


Combining Proposition 5.3 with Theorem 5.1 we have established the claimed non¬ 


degenerate integrated decay with derivative loss result for the wave equation. 


5.2. Proof of Theorem |1.1| for Spin 1. The proof of the theorem for the Maxwell field 
follows a similar pattern to that of the conformal scalar field. There is a simplification 
owing to the fact that the energy-momentum tensor is trace-free, and the elliptic estimate 
takes a slightly different form. 


Proposition 5.4 (Boundedness of energy). Let T 2 > Ti. Suppose that F is a solution 
of Maxwell’s equations 29), subject to the dissipative boundary conditions (fO) as in 


Theorem 4-2. Then we have: 



V I if 1 r'^drduF 

" ■^00 

-|- |ii|^^ r'^drduj. 


[Ti,T 2] 


(^\r'^EAf + |r^iiA|^^ dtduj 
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Proof. We apply Lemma 2.1 to the vector field = '^[F]abT^- Consider the term on 

J^. We have 


, fr r 1 = lim r^(l + EaH sdidw 

r—>-oo 1^01,T2] r—>-oo /fi™iT2l 

t/ tJ Z-if' 


lim 


'^00 

t'^HaI^ dtdio 

dtdio. 


= — lim 


r^oo /^[T'iT2l 


2 hlTl,T2 


Here we have used the dissipative boundary conditions (40). Since Va^JI“ = 0, there is no 
bulk term, and a simple calculation gives 


^JIan“dS'Er = 
/Et 


1 


“''^2 9 , 


+ r’^drdio, 


which completes the proof. 


□ 


We next show an integrated decay estimate with a loss in the weight at inhnity: 

Proposition 5.5 (Integrated decay estimate with loss). Let T 2 > Ti. Suppose that F is a 
solution of Maxwell’s equations 29), subject to the dissipative boundary conditions (fO) 


as in Theorem 4-2. Then we have: 


‘^ + drdujdt + 

" Wl + r^ 

<3 j r'^drdoj. 


lf,lTl,T2] 

■^00 


r‘^Ef + ) dtduj 


Proof. We apply the divergence theorem to integrate the current "’^1“ = over <S'[Ti,r 2 ]' 

Recalling the expression 0 for ^TT, and noting the fact that the energy-momentum tensor 
is traceless, we have 


Div = ^TTabT^^ 
Too 


VTw 

\E\'^ + \H\- 

vT+T^ 


RT^Vl + r-^ 
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Now consider the flux through a spacelike surface We have 


X TT 


= / Tor 


/St 


Vl + T- 


-.drdu! 


= I e^^EAHB 




-drduj 


so that if t G [Ti,r 2 ] we have by Proposition 5.4 


'St 


1 

< - 
“ 2 


f (\E\‘^+ \H\Ar‘^drduj <- f (lEl"^ + r^drdu 

JTit ^ 


Next consider the flux through a surface of constant r. We have: 


L 


[Ti,T 2\ 


^ Ia'fn°‘dS^ = 


jfp^l,T2\ 


\/l + r'^dtdu 


[ (- \Er\^ - \Hr\‘^ + \Ea\^ + \Ha\^) rWl+r^dtduj 


so that 


lim / ^Sa'm^'dSy < — / 

1 —^■OO lylTi,T2] Sr |■y[Tl,T2] 

O 2 -ir 'J ^oo 


+ 2 


/St]^ 


dtdu 

(\E\^ + \H\^^ r'^drdoj, 


where we use Proposition 5.4 to control the angular components of E, H at infinity. In¬ 
tegrating Div over <S'['r^^r 2 ]! ^PPly^g the divergence theorem and using the estimates 
above for the fluxes completes the proof of the Proposition. □ 

We next seek to improve the weight at infinity at the expense of a derivative loss. The 
first stage in doing this is an elliptic estimate. First note that by commuting the equations 
with the Killing held T we have: 

(47) f {\dtE\^ + \dtH\^) -jL=drdLodt + {\r^dtEf + \r^dtHf) dtdu 


(48) 


'%lT2l 

< 3 


(^{dtEl"^ + \dtH\^'^ r'^drdoj. 


We can use the evolution equations to control the right hand side in terms of spatial 
derivatives in the slice if we choose. We will require the following Lemma: 


Lemma 5.5. Let K he the veetor field 
2r 


K := 


\/rT 


-.Hr 


{rf^HA) er-^ (rVl + ca 


1 + r- 


^H^er- 
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Then if H satisfies the constraint equation (29) we have the identity 
2 


\/rT^ 


zdr {ry /1 + r'fj^ 


Hr 


= 2 






+ 


(1 + r^) 2 


\Hr\ +Div K 


Proof. See Appendix 7.2 


□ 


Now consider the Maxwell equation (25): 
r 


VTT' 


-,dtEA = eA 


B 


dr (r\/1 + r’^Hs) — r'^ 


squaring this and multiplying by (1 + r^) 2 ^ we have 


\dtEA\^ 

(1 + r2)2 


1 




r^AHrV + 


a 

dr 


Vl + t^Ha) 


(49) - dr (r\/ 1 + r‘^HA \ r'^^ 

V1 + V / 

From here we readily conclude: 


Hr. 


Theorem 5.2 (Higher order estimates). Suppose that E is a solution of Maxwell’s equa¬ 
tions, as in Theorem 14.21 Then there exists a constant C > Q, independent of Ti and T 2 
such that we have 


1 






+ 2 


\/rT 


r-^AHrV + 




a 

dr 


A 

[rVl + t^Ha) +\rfAHBf 


+ H E yr^dtdrdoj 


< C 


' Stj 


+ \H\^ + \dtE\^ + \dtH\^^ r^drduj. 


Proof. We integrate the identity (49) over 5'[ri,T2]- We control the left hand side by the 
time commuted integrated decay estimate (47) since r^(l + r^)“^ < 1. The right hand 
side, after making use of Lemma 5.5 and applying the divergence theorem will give us 
good derivative terms, a good zero’th order term which we can ignore and a surface term. 


The surface term at infinity give: 


16 


a term proportional to , integrated over the 


cylinder, which we control with the estimate in Theorem |5.5[ This immediately gives the 


^®Notice that there is no surface term contribution at infinity arising from ( P due to 

■\J l+r^ 

the decay of the angular derivative of iF in r by the local well-posedness theorem (Theorem | 4 . 2 |. 
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result for all of the terms except the term , which we obtain from a standard 

elliptic estimate on the sphere (see for instance Proposition 2.2.1 in [1]) : 


\fAHB\^<C [ 

JS2 


AB 


^aHbI^ + \f Ha\^ 


after noticing that we already control aHb and ^ Ha with a suitable weight from 
(24) and (29). Finally, we note that the estimate can be derived in an identical manner for 
E. □ 


Theorem 5.3 (Full integrated decay). Suppose that F is a solution of Maxwell’s equations, 
as in Theorem\4.^ Then we have 


i 

\Er\^ + \Hr\^ + \Ea? + \Ha\^ 

'^[Ti,T2\ 



< C 


-Ti 


r'^dtdrdoj 

+ 


Efr 

2 

+ 

Ea 

2 

-h 

Hr 

2 

-h 

Ha 


r^drduj 


for some C > 0 independent ofTi and T 2 . 


Proof. We combine the result of Theorem 5.2 with the Hardy estimates of Lemma 5.4 with 
a = 1 to improve the weights near inhnity in the integrated decay estimates, making use 
of a cut-off in much the same way as for the spin 0 problem. □ 


5.3. Proof of Theorem |1.1| for Spin 2. As a useful piece of notation we first introduce 
the trace-free part of Eab by 


Eab — Eab — -^^abEc^ — Eab + -^^^abEt 


and similarly for Hab ■ The boundary conditions of Theorem 
expressed as 


4.3 


may then be conveniently 


r'^ [Eab + ^{A^Hb)c 


0 , 


as r —00 , 


and we also have \Eab\‘^ = \Eab\‘^ + \ \Eff\^■ We will first prove the boundedness and the 
degenerate integrated decay statement of the main theoremj^ 

Proposition 5.6 (Boundedness of energy and integrated decay estimate with loss). Let 
T 2 > Ti. Suppose Eab and Hab are solutions to the spin 2 equations subject to (41) as in 


17 


As we will see in the proof, the key difference to the case of the wave- and Maxwell’s equations is that 
here we will have to establish both these estimates at the same time! 
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Theorem \4-3[ there exists a constant C > 0, independent ofTi and T 2 such that we have 
{EabE^^ + rWl + r^dtdrduj + j 1 - 2(1 + r‘^)drdu: 

+ {EabE^^ + r^dtdiw 

J Soo 

<C [ {EabE^^^ + HabH^^} r2(l + r^)drd 0 J. 

JEtj '' 

Proof. Let us now introduce the vector field 


d S[Ti,T2\ 
(50) 


1 


Y = T+^X. 
\/3 


We certainly have that Y is timelike, since 


g{Y, Y) = -(1 + r2) + = -1 - < 0. 


Moreover, we have 


Y 


vr = 


%/3 




The current that we shall integrate over is 

J“ = QIJP^Y^Y^. 


Clearly, we have by Lemma 3.1 that the flux of J through a spacelike surface with respect to 
the future directed normal will be positive. Moreover, since Q is trace-free and divergence 
free, Div J will also be positive. To establish a combined energy and integrated decay 
estimate, we simply have to verify that the surface term on ^ has a definite sign (and 
check the weights appearing in the various integrals). We shall require some components 
of Q, which are summarised in the following Lemma: 

Lemma 5.6 (Components of Qabcd)- With respect to the orthonormal basis in which we 
work, we have 

Qoooo = \Eab\^ + \Hab\^ + 2 + 2 \HAr\^ + iL-rrl^ + \Hrr\^ 

Qooor = 2 {EAce^^HB^ + EAfe^^HBf) 

QoOfr = \Eab\^ + - \Erf\^ - \Hrf\^ 

Qorrr = 2 {EACe^^H b^ - EAr^^^HBr) 

Qrrrf = IEaB^ + “ 2 \EArf “ 2 \HArf + \Errf + \Hrrf 

With our definition of the current J above, we apply the divergence theorem. Lemma 


2.1 We now verify that all the terms have a definite sign. 
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a) Fluxes through We compute 


/ = / Q(eo,y,y,y) 

'T.t JT.t 


VT+~r 


zdrduj. 


Now, defining F = (1 + 4r^) ^Y, we have that 


-gieo,y) = 


I 1 + r^ 

1 + |r2 


< \ - 


so by Lemma 3.1 we deduc^ 


so that 
and hence 


Q(eo,F,F,F)~goooo, 
Q{eo,Y,Y,Y)r.{l + r^yQoooo, 
JarfdS^, ~ ^ (\Eab\^ + r2(l + r‘^)dr. 


b) Bulk term We have 

DivJ = = 




\/l + r' 


:Q(eo,eo,F,F). 


Again applying Lemma 3.1 to F and rescaling, we have 

DivJ ~ (1 + r^) 2 Qqqoq. 

As a result, we have 




Div Jdr] ~ [ (^\Eab\‘^ + r^-\/l + dtdrdoj. 


’[ri.T2l 


c) Boundary term at Finally, we consider the flux through surfaces We have 

f Ja'm°'dSf^^= f Q{er, F, F, Y)r‘^\/l + r'^dtduj, 

Now, inserting our expression for F and expanding, we have 
Q{er,Y,Y,Y) 


— (l + 2 ^g^rgoo + "v/s 

so that 


7* 1 7*'^ 

Qt^OO “1“ I nQrrrO “1“ “ 7= ^Qrrrr 


\/l + r' 


1 +r^ 


3\/3(i + r2)2 y 

lim Qipri F, F, F)r 1 A = lim r | grOOO F V^grrOO F QrrrO F '^Qrrrr ) • 
•^oo r^oo 3.^3 y 


^®We write / ~ g to mean that there exists a numerical constant C > 0 such that C ^ f < g < Cf 
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Let us consider the first and third terms on the right hand sid^^ We have: 
QrOOO + QrrrO = Hb^ 

= 2 (^Eab - \sabEc^ + €(a^^b)c) (e^^^ - 
-2\Eab? -2\Hab? 


so that 


lim r® (Qrooo + Qrwo) = -2 lim r® {\Eab\^ + \Hab\‘ 

r^oo r^oo \ 


where we make use of the boundary condition. Taking this together with the expressions 
for f^T’T’OO? TTTT v^e have. 


lim Q{er, Y, Y, Y)r‘^\/l + = — lim r® 





so that 


and 


lim Q(e^,y,F,F)rVl + r2~- lim + \Hab\‘') 

r—^oo r—^oo V / 


■■,[Ti,T2\ 


Ja'm'^dSf. ~ — 


ff,[T’i,T2l r^oo 

^OO 


li^r® (^\Eab\‘^ + dtduj. 


Taking all of this together, we arrive at the result. 


□ 


As we did in the Maxwell case, we shall now use the structure of the equations to 
allow us to establish (weighted) integrated decay estimates for all derivatives of the fields 
E, H. To control time derivatives we can simply commute with the Killing field T and 
apply Proposition |5.6[ To control spatial derivatives we replace the time derivatives by the 
equations of motion and integrate the resulting cross terms by parts making use also of the 
constraints equations. The remarkable fact is that in the process we only see spacetime 
terms with good signs and lower order surface terms that we already control by the estimate 
before commutation. 

We will note the following useful result, which allows the cross term to be integrated by 
parts: 

Lemma 5.7. Let K he the vector field 

(51) K:= (2r^f^HBcH^r- , {Hb^A Cr - [r{l + r^)H^^] eg. 

V VI + J yjl + r^ or ^ 


factor ^ appearing in the definition of Y was chosen to arrange a cancellation between these 


terms. 














ASYMPTOTIC PROPERTIES OF LINEAR FIELD EQUATIONS IN ANTI-DE SITTER SPACE 


34 




If H satisfies the constraint equation (Con Hb), we have the identity 
2r^ ,c 

DivK = - 

Proof. See Appendix |7.3[ □ 

Proposition 5.7 (Higher order estimates). Let T 2 > Ti. If Eab,IIab solve the spin 2 
equations subject to the dissipative boundary condition (41) as in Theorem 4-3. then we 
have 


^[Ti,T2] 


|9r(T(l + r‘^)IlBc)^ + lltZL {t'^IIbt)^ ^drdtdoj 


1 + r 


+ 


\. 


\ 


^[Ti,T2] 


r^\1lBHcr? + r^\1/AHBC? + (H o E)|drdtda; 


<C f {EabE'^'^ + HabH^^ + EabE^^ + Habll^^^] r"(l + r^)drdiw, 

j Stj 


for some C > 0 independent ofTi and T 2 . 

1 d 


Proof. Recall now ( |Evor EabV - 
r OEab 


vTT^ dt 


= e 


{A 


c 


Consider 




Vi + 

1 d 


(r(l + r‘^)HB)c) - 'r"^\c\HB)-T 


[Vm^ar (W + r^)HBc)-rfcIiB, 

' ® (r(l + r‘)H‘>B) - rf’^HBf 


\/l + r2 dr 
Vl + 1-2 dr 


Hrr ] - rf^HB 


drir^Hrf) - rf^HBr + 

2 Hrr 


VTTr^ 


Vl + 

Now, since for any 2—tensor Z on we have Z]^ab] = Zcd), we deduce that if 

the constraints hold then (EvoP Eab) may be re-written: 

r OEab 


\/l + r^ dt 
whence we deduce 

Xab := ^ ^ 


= ^A 


C 


1 d 
Vl + dr 


(r(l + r‘^)HBc) - rfcHBf 


^AB 
\/l + r' 


-Hrr^ 


Vl + r2 dr 


(r(l + r^)HBc) - rfcHBf = e^c 


VT+~r2 dt Vl + 
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Now, using the second equality and applying the estimates in Proposition |5.6| for H and 
the commuted quantity E, we can verify that 


Xdtdrdoj 


'^[Ti,T2] 


< c 


< c 


hTi.Tj] 


Eab 


'St, 


[EabE^'^ + + E, 


+ r^\/l + r'^dtdrdio 

bE^^ + HabH'^^} r^{l + r‘^)drdu. 


We also have, however, 

f X abX^^ dtdrduj 

''^[Ti,T2\ 


'^[T-l,T2\ 


1 _j_ 


\dr{r{l + r^)HBc)\^+ r^\1lBHa 


[r(l + r^)H^c] \dv 


\/l + r 


'^ITi,T2] , 

+ 2 EE 


1 _j_ 


\dr{r{l + r^)HBc)\^+ r^\fBHa 


\dr{r^HBr)f \drj- \r^HBrf dtdu, 

I Lioo 


where in the last step, we have used the result of Lemma 5.7 to replace the cross term with 
a good derivative term and a surface term. 

It remains to control the term aHbc\'^■ Notice that by (EvoL Eat) and (Con Hb), 
we have 


Ae^^1/BHcA\^+ A1/ Hab? <c[ -^\EAr? + 

\ 1 + 


1 + 


d 




+ 


\HA-r\^ \ 

1 _|_ ^2 I • 


On the other hand, since 

e^oe^^^BHcA = f^HAD - fnHA^, 


we can apply the standard elliptic estimate (see for instance Lemma 2.2.2 in ID) 

I <C f (\e^^fBHcA\^ + If'^HABl^ + IfAHB^l^) 

JJ 

to obtain the desired bounds for AdlBc\^- 

This gives all the desired estimates for the derivatives of 77. As in the Maxwell case, 
similar bounds for the derivatives of E can be derived in an identical manner. □ 
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Finally, much as in the Maxwell case, we apply the Hardy inequalities to establish 
integrated decay of the non-degenerate energy with the loss of a derivative: 

Theorem 5.4 (Full integrated decay). Suppose that W is Weyl tensor, satisfying the 


Bianehi equations with dissipative boundary eonditions, as in Theorem f.3. Then there 
exists a constant C > 0, independent of T such that we have 


'S[Ti,T2] 


EabE^^ + 




<c { 


EabE^^ + + EabE^^ + HabH^^ } ^’^(1 + r^)drduj. 


Proof. Using the result of Proposition 5.6, 5.7 with the Hardy estimates of Lemma 5.4 


to improve the weights near inhnity in the integrated decay estimates, making use of a 
cut-off in much the same way as for the spin 0 and spin 1 problems, we obtain the desired 
estimates for \Eab\, \Hab\, \EAf\ and \HAf\- Finally, the bounds for \Erf\ and |Hff| are 
obtained trivially using the trace-free condition for E and H. □ 


5.4. Proof of Corollary 1.2 (uniform decay). In this subsection, we show a uniform 


decay rate for the solutions to the confomal wave. Maxwell and Bianehi equations, hence 
proving Corollary |1.2[ We will in fact prove the uniform decay estimates for all of the 
equations at once by showing that this is a consequence of the bounds that we have obtained 
previously. More precisely, we have the following general lemma: 

Lemma 5.8. Let 'L 6e a solution of either the conformal wave, Maxwell or Bianehi equa¬ 
tions. Suppose that we are given some positive quantity <S['I'](t) depending smoothly on 'h 
and its derivatives at some time t which satisfies: 

1. <f[4'](t) is a non-increasing function oft, 

2. For every 0 <Ti < T 2 , <f['L](t) satisfies the integrated decay estimate: 


rT2 


/Ti 


£[m)dt <C{£[mTi)+S[dtmT,)} , 


for some C > 0 independent ofTi and T 2 . 
Then we have the estimate 


£:[4/](t)< 


c„ 


{1 + ty 


( 0 ) 


fc =0 


for some constants C„ > 0 depending only on n and C. 
Proof. Let us set 

£^’‘\t)=£\{dt)’^'^] {t). 













ASYMPTOTIC PROPERTIES OF LINEAR FIELD EQUATIONS IN ANTI-DE SITTER SPACE 


37 


We calculate 

(1 + t - + y* ^ ((s - ds 

= S^°\t)+ [ S^^\s) + {s - Ti)S^^\s)ds 

Jti 

<£^^\t)+ [ £^^\s)ds, 

Jti 

where we have used the monotonicity of £^^^ to obtain the last inequality. Now, it follows 
from the assumptions of the Lemma that 

£^^\t) + j £^^\s)ds<Ci[£^^\Ti)+£^^\Ti)^ , 

which together with the preceding estimate immediately imply 

(52) f(0)(t) < . 

Taking Ti = 0, this in particular implies the conclusion of the Lemma in the case n = 1. 

To proceed, we induct on n. The n = 1 case has just been established. Suppose now 
that the statement holds for some n. Noticing that the equation commutes with dt, we use 
the induction hypothesis for both T and dt^ to obtain 




a 


n+l 


(i + t)' 


^£:W(o). 


fc =0 


Now, we apply (52) with Ti = | to deduce 


n+l 


< 


2 

ClCn 


< 


(1 + |)-+1 

ClCn2 
(1 + 


2 ^ fc=0 
,n ^+1 




k=0 


whence the result follows. 


□ 


Proof of Corollary 1.2. For the conformal wave and Maxwell equations, Corollary 1.2 fol¬ 
lows immediately by applying Lemma |5.8| to the quantity 


mit) = 


/ , r^drdu) 

hr Vl + r^ 


which is monotone decreasing and satisfies an integrated decay statement with loss of one 
derivative. For the Bianchi equations, this quantity is not monotone decreasing (merely 
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bounded by a constant times its initial value). We can circumvent this by applying Lemma 

S[W]{t)= [ Q{eo,Y,Y,Y)- 


5.8 to the quantity: 


'St 


V^' 


-drduj, 

i > 


which is monotone decreasing and satisfies an integrated decay statement with loss of one 
derivative. Noting that 

S[W]{t) ~ f f r^drduj, 

Jut V 1 + 

we are done. □ 


5.5. Proof of Theorem |1.3| Gaussian beams. It is noteworthy that in the first in¬ 
stance, for all of the integrated decay estimates we obtained above the r-weight near infinity 
is weaker than that for the energy estimate. In particular, in order to show a uniform-in- 
time decay estimate, we needed to lose a derivative. In this section, we show that without 
any loss, there cannot be any uniform decay statements for the conformal wave equation. 
Moreover, an integrated decay estimate with no degeneration in the r-weight does not hold. 

In order to show this, we will construct approximate solutions to the conformally coupled 
wave equation for a time interval [0, T] with an arbitrarily small loss in energy. We will in 
fact first construct a Gaussian beam solution on the Einstein cylinder and make use of the 
fact that (one half of) the Einstein cylinder is conformally equivalent to the AdS spacetime 
to obtain an approximate solution to the conformally coupled wave equation on AdS. 

In the following, we will first study the null geodesics on the Einstein cylinder. We 
then construct Gaussian beam approximate solutions to the wave equation on the Einstein 
cylinder. Such construction is standard and in particular we follow closely Sbierski’s geo¬ 
metric approach m (see also 1271 EH]). After that we return to the AdS case and build 
solutions that have an arbitrarily small loss in energy. 

5.5.1. Geodesics in the Einstein cylinder. We consider the spacetime {Me, 9e), where M. e 
is diffeomorphic to M x and the metric is given by 

(53) qe = —dt^ + d'lp'^ -|- sin^ •il){d9‘^ + sin^ Odcj)'^). 

We will slightly abuse notation and denote the subsets of Me with notations similar to 
that for the AdS spacetime. More precisely, we will take 

YT-.= {{t,^,e,cP)-.t = T}, 

^ITi,T 2] ■= {{t, - TiYt < T 2 ]. 

Take null geodesics 7 : (— 00 , 00 ) —)■ AIe in the equatorial plane {9 = |}. In coordinates, 
we express 7 as 

{t,'ip,9,(l)) = (r(s),T(s), ^,$(s)). 

Since ^ and ^ are Killing vector fields, E and L defined as 

t = E, sin^ Ti = L 
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are both conserved quantities. Here, and below, we use the convention that ’ denotes a 
derivative in s. We require from now on that 0 < |L| < E. The geodesic equation therefore 
reduces to the ODE 


(54) 


^ - 


L 2 


sim T 


Solving (54) with the condition that T achieves its minimum at s = 0 and r(0) = 4>(0) = 0, 
we have 


T{s) = Es, 

. , ^{E^-L^)sm\Es)+L^ 

smT(s) = 4^4-^^ 

and 

a./ ^ _ r 

Jo (E2-L2)sin2(Es')+ ■ 

In order to invert the sine function to recover T, we will use the convention that for 
Es G [2kTr — |, 2k7r + |), and A: G Z, we require 'h(s) G [0, |]; while for Es G [{2k + l)7r — 
(2k + l')7r + j), and A; G Z, we require 'h(s) G [f ,71]. Notice that this choice of the inverse 
of the sine function gives rise to a smooth null geodesic. 

Moreover, direct computations show that 


t(s) = E, 4/(s) 


EVE 2 - L 2 sin(Es) 
x/(E2 - L2)sin2(Es) + L2’ 


4>(s) 


LE^ 

(E2 -L2)sin2(Es) +L2' 


5.5.2. Constructing the Gaussian beam. Given a null geodesic 7 on {AiEidE) as above, 
we follow the construction in Sbierski HZl to obtain an approximate solution to the wave 
equation on {Me^Qe) which is localised near 7 and has energy close to that of 7 . We 
first define the phase function y? and its first and second partial derivatives on 7 and then 
construct the function y? in a neighbourhood of 7 . We also define the amplitude a on 7 . 
More precisely, on 7 , we require the following conditions: 

( 1 ) <y9(7(s)) = 0 

( 2 ) dip{j{s)) = 7 ( 5 ) 1 , 

(3) The matrix M^i, := d^di,g:i{'^{s)) is a symmetric matrix satisfying the ODE 

= -A- BM - MB'^ - MCM, 
as 

where A, B, C are matrices given by 

= ^{d^dpg^'')d^ipd^(p, 

Bnp = d^gff^dpcp, 

C.p = g^P, 

and obeying the initial conditions 

(a) M(0) is symmetric; 

(b) M{Q)p,r = {dp^m-, 
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(c) ( 8 ) dx'^ is positive definite on a three dimensional sub¬ 

space of T^(o)-^ that is transversal to 7 . 

(4) 0 ( 7 ( 0 )) / 0 and a satisfies the ODE 

2 grad(/?(o) -|- □(/? ■ a = 0 


along 7 . 

The results in m ensure that , d^ip\^ , d^dy(p\^ , a\^ can be constructed satisfying 
these conditions. We then let to be a smooth extension of (p away from 7 compatible with 
these derivatives. Likewise a/p is defined to be an extension of a\^ as constructed above. 
Moreover, we require aj\f to be compactly supported in a (small) tubular neighborhood M 
of the null geodesic 7 . 

We define the energy for a function on by 


EtH 


1 

2 


{{dtwf + {d^wf + 


(pewf^ 

sin^ -ip 


H--—^^ 7 ^) sin^ ipd'ipduj. 

sin Ip sin 6 


Here, as elsewhere, doj = siwOdOdcp is the volume form of the round unit sphere. We also 
associate the geodesic 7 with a conserved energy 


E{l) = -gEii,dt). 

Notice that this agrees with the convention E = T used in the previous subsection. 

The main resullj^ of Sbierski regarding the approximate solution constructed above is 
the following theoreirj^ (see Theorems 2.1 and 2.36 in [IZ!): 

Theorem 5.5. Given a geodesic 7 parametrized by E and L as above, let 


where ajp and (p are defined as above. Then obeys the following conditions: 

(1) \\E\we,l,xm\\l2(S[o,t]) - 

(2) Eo{we,l,x,jv) 00 as a 00 ; 

(3) we,l,x,jV is supported in M, a tubular neighborhood of^; 

(4) Fix p, > 0 and normalize the initial energy of we,l,x,jV by 


WE,L,xpf 


WE,L,X,Ar 


Eo{we,L,X,PJ') 


•E(7). 


^^Translated into our notation, the result in m requires the following bounds on the geometry of 
{ME,gE)- 

-C < g{dt,dt) < c < 0, \Vdt{dt,dt)\ + \Vdt{dt, eOl -t \Vdt{ei, efi\ < C, 
where Ci is an orthonormal frame on the slice. These estimates are obviously satisfied in our setting. 

^^Regarding point (4) in the theorem below, the original work of Sbierski gives a more general charac¬ 
terization of the energy of Gaussian beams in terms of the the energy of geodesics on general Lorentzian 
manifold. Since in our special setting, the energy of a geodesic is conserved, we will not record the most 
general result but will refer the readers to nzi for details. 
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Then for M a sufficiently small neighborhood of 7 and A sufficiently large, the 
following bound holds: 


sup 

te[o,T] 


Et{wE,L,XM) - ^( 7 ) 


< 


We also need another fact regarding the second derivatives of ip which is a consequence 
of the construction in |T7j (see (2.14) in [T7]l: 


Lemma 5.9. |^) = ^{Vip |^) = 0. Moreover, 9‘(VV(/? |^) is positive definite on a 

2>-dimensional subspace transversal to 7 . 

The fact that app is independent of A together with Lemma |5.9| imply that the Gaussian 
beam approximate solution constructed above has bounded norm independent of A. We 
record this bound in the following lemma: 


Lemma 5.10. Let we,l,\jV be as in Theorem 5.5. 


The following bound holds: 


We,L,\M\l'^(S[o^t\) - ^(^)- 


5.5.3. The conformal transformation. Once we have constructed the Gaussian beam for 
the wave equation on the Einstein cylinder, it is rather straightforward to construct the 
necessary sequence of functions using the conformal invariance of the operator 

L = Og- ^R{g), 

where R{g) is the scalar curvature of the metric g. 

We first set up some notations. We will be considering M.e restricted to i/; < | as a 
manifold with boundary diffeomorphic to M x S^. The interior of this manifold will also be 
identified with M.AdS via identifying the coordinate functions {t,6,4>), as well as 


tanijj = r. 


It is easy to see that gE and gAdS are conformal. More precisely, gE as before can be 
written as 

Pe = —df^ + df'^ + sin^ fi{d6‘^ + sin^ Odffi)-, 
while in the {t, 9, fi) coordinate system, pAdS takes the following form 

9AdS = —Vt ( — dt"^ + dfi'^ + sin^ if{d9‘^ + sin^ Odcjfi)). 
cos^ yj 

We now proceed to the construction of the approximate solution to the conformally 
coupled wave equation on AdS using the conformal invariance of L. More precisely, we 
have 


Lemma 5.11. Let w be a function on Me restricted to tp < 5. Then we have 


1 


(1 + r^) 2 


{Og^W -W)= Dg^^s 


w 


\/r 


+ 2 


w 


VT+~r2 


Proof. This can be verified by an explicit computation. 


□ 
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Figure 1. An typical null geodesic 7 shown in black projected to a surface 
of constant 9 in the optical geometry of the Einstein cylinder (left). The 
corresponding curve for the anti-de Sitter spacetime, after reflections at 
(right). 


We are now ready to construct the approximate solution. To heuristically explain the 
construction, consider Figure On the left we sketch the curve 7 along which our ap¬ 
proximate solution is concentrated in the Einstein cylinder. To allow this to be plotted, 
we project onto a surface of constant t, 9 which carries a spherical geometry and can be 
visualised via its embedding in M^. We actually wish to construct an approximate solution 
on anti-de Sitter, which is conformal to one half of the Einstein cylinder. To do so, we 
restrict our attention to one hemisphere and arrange that whenever the curve 7 strikes the 
equator of the Einstein cylinder, it is reflected. Each time this occurs, we shall arrange 
that the Gaussian beam is attenuated by a factor which depends on the angle of incidence 
(which in turn depends on E and L). The more shallow the reflection, the weaker the 
attenuation. Crucially, the time, t, between reflections is independent of the angle that 
we choose. As a result, by taking the angle of incidence to be sufficiently small, we can 
arrange that an arbitrarily small fraction of the initial energy is lost at the boundary over 
any given time interval. 

To be more precise, let us fix T > 0. We will construct an approximate solution for 
t G [0, T]. Take iV G N be the smallest integer such that T < GIII+ITe, wig then construct an 

approximate solution for t G ( —f, iM+TzEj starting from the function we,l,x,M constructed 
previously. This will correspond to a Gaussian beam which strikes the boundary ~ 2N 
times between t = 0 and t = T. Notice that the geodesic 7 has the property that it lies in 
the hemisphere {^p < for t G (2/^71 — |, 2A:7r-|- |) (for k G X) and that it lies in the other 
hemisphere, i.e., {i/j > |}, for t G {{2k + l)vr — {2k + l) 7 r-|- |) (for k G Z). Therefore, we 

will assume without loss of generality that the neighborhood M has been taken sufficiently 






ASYMPTOTIC PROPERTIES OF LINEAR FIELD EQUATIONS IN ANTI-DE SITTER SPACE 


43 


t 



i, = ,r/2 


Figure 2. The cutoff functions. 


small such that it lies entirely in {il) < for t G (2/c7r — |^,2fc7r + (and entirely in 

{ip > |} for t G {{2k + l)7r — {2k + l)7r + |)), where k G Z. 

Now for {t, Ip, 9, p) \ii M.E restricted to we define for t G [0, T] 


(55) 


N 


UE,L,\{t^'P^&}P) ^ X;:g(2fc7r-^,2fc7r+^] 


fc =0 

Af 


Vl + 




7r-^,(2fc+lK+^]- 


fc =0 


WE,L,X,jv{'t, TT - V', 6*, <)>) l{V)<f} 
VTTr^ 


where i? is taken to be i? = — and x is the indicator function. Notice in 

particular that when the time cutoff function is 0 in the first term (resp. in the second 
term), the support of we,l,\jV is entirely in {ip > |} (resp. {ip < |}). We also depict this 
in Figure where we denote 

N N 

Xl{t) = ^Xte(2fc7r-^,2fc7r+^]> X2{t) = Xigf(2fc+lW-gf .f2fc+lW+gf 1' 

k=0 k=0 

The definition above is such that in the time interval t G (—X’ xl’ we,l,\jG 

constructed previously, restrict it to < | and rescale it by Then on the time 

interval t G (tt — ^,7r + we take the part of we,l,\jG that is supported in > |, 
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reflect it across the lA = f hypersurface, rescale by a factor and then multiply by the 

factor R. As we will see later, the factor R is chosen so that the boundary conditions are 
approximately satisfied. We then continue this successively, taking parts of the solutions 
in lA < f and ip > reflecting when appropriate, and multiplying by factors of i?’s. 


Lemma 5.12. 

(56) 


The function ue,l,x defined as in (55) has the following properties: 
Eo[ue,l,\\ ^ oo as X 


oo; 


(57) 

(58) 


{Eig^as^E,L,\ + 2uE,L,\fr‘^ + r^) dr du) dt < C{T)] 


^[0,T] 


[_,^^Xdt{ruE,L,\) + r^dr{ruE,L,\))‘^ du dt < C{T). 


Proof. First note that Eq[we,l,x] < CEq[ue,l,x], so the first claim follows from Theorem 
5.5 Next recall that we have assumed that the neighbourhood AA has been taken sufficiently 


small so that for each of the summands, the indicator fnnction y is constant on the support 
of we,l,x,jV when restricted to a hemisphere. It therefore suffices to consider only the 
contributions from we,l,x,jV (as no derivatives fall on the indicator functions). We will 
write 


ue,l,x = ue,l,x,i + ue,l,x,2, 


where ue,l,x,i is the first snm in (55) and ue,l,x,2 is the second sum in (55). For ue,l,x,i, 
we apply Lemma 5.11| to get 


< 


*[0,T] 

I. 


{E\g^^gUE,LXi + 2uE,L,x,i)‘^r‘^{'i- + r^)drdLjdt 

{E\gj^WE,L,x,jV - WE,L,x,jv)^^^i'^ + ^^) dr duj dt 




{Ugj^WE,L,x,N - WE,L,xjsr) sin^ Ip dip duj dt 


5[0,T] 


<C{T). 


In the last line, we have used the estimates in Theorem 5.5 and Lemma 5.10 By a 
straightforward identification of the two hemispheres in AI^;, we can prove a similar bound 
for ue,l,x,2- 

/ (□ 5 Ads^i?T,A ,2 + 2 uE,L,x, 2 )‘^r‘^{l + r^) dr du dt 

< / {ngj^WE,L,xM-'<^E,L,xMf^^^‘^'^dipdudt<C{T). 

Sin 'n 


This concludes the proof of the second claim. 

We now show that the bonndary terms are appropriately bounded in L^. By construc¬ 
tion, there are only contributions to the boundary terms in a neighborhood of :^(t — ^) G N. 
Since there are at most 0{N) = 0{T) such contributions, it suffices to show that one of 
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them is bounded. We will look at the boundary contribution near t = f, which takes the 
form 

lim {rdtUE,L,\ + r‘^dr{ruE,L,\)){t, r, 6, (f) 

r^oo 

= lim {dtWE,L,\^ + d,ijWE,L,\^){t,'4’,0,(t)) 

(59) + lim R{dtWE,L,\,M + d'>pWE,L,\j^){t,'n - 

=i\aM{{dt^ + + R{dtv - 

+ {{dtaM + d^a-M) + R{dtajv - 


The latter term is clearly bounded pointwise independent of A. The first term has a factor 
of A and we will show that it is nevertheless bounded in since by the choice of R, 
{{dt^p + d^(p) + R{dtp> — d^if)) vanishes on 7 . More precisely, by points (1), (2) in Section 
|5.5.2| and Lemma |5 .91 we have 




Ett 

2 


LE^ 


lo {E^ -L‘^)sm\Es') + L‘^ 


ds 


/\ 2 \ 


for some a > 0. We further claim that 


+ d^ip) + R{dtp> - d^(p)) 


{t 2>b 2 T 2’*^ fo^ (E'-i-L'^)sin^(E3’)+L'^ ^ 

This follows from the fact that d(p = 7 t, on 7 and the choice of R. More precisely, we have 

i{dt^ + d^ip) + R{dt^-d^ip)) I _ ^ ^^2 , 

it-2-V-2.'^-2’'P-Jo (E2-L^)3in2(Es')+L'2> 

= {-E + ^E^ - L2 + R{-E - y/W^I?)) 

= 0 . 

The desired bound on the boundary then follows from Lemma 


= 0 . 


5.13 


□ 


Lemma 5.13. Suppose that f is a function defined on {ip = |} which vanishes to order 0 
at{t=^,e = ^,cP = fo~ (]^2ZI2)S(eT)W)- 

f dco dt\fe‘^^'^\‘^ < 


Proof. In order to simplify notation, we define 

LE"^ 

■“ Jo {E^ -L‘^)siifi{EP) + L^^^ ' 
The statement that / vanishes to order 0 is equivalent to 

I/I < C (^{t - -)^ + {6 - -)^ + (</>- 4>o)‘^^ 
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Thus 


dujdtlfe^^'^f = J dudt\ff 


r*27r 






<C I sinOdO I dcj) I dt l{t —+ {9 —+ {(j) — e 

lo Jo Jo ^ 2 / 

,|2 . 5 


_ ^ n 2\ 2Aa((t-f)2 + (0-f)2 + (0-</,o)2) 


<c 


\x 


2 g- 2 aA|a,| < ^^- 5 ^ 


where in the last line we simply scale A out of the integral. 


□ 


5.5.4. Building a true solution. Given the approximate solution constructed above, we are 
now ready to build a true solution to the homogeneous conformally coupled wave equation 


with dissipative boundary condition. To this end, we need a strengthening of Theorem 5.1 
which includes inhomogeneous terms: 

Theorem 5.6. Letu be a solution of the inhomogeneous conformally coupled wave equation 
(60) + 2u = f in AdS 

with finite (renormalized) energy initial data 

Et^ (m) < oo 

and subject to inhomogeneous dissipative boundary conditions 


d{ru) 2'9(tu) 
dt dr 

Then we have for any Ti < t < T 2 : 


as r —>■ 00 . 


Et[u] < Cti,T2 Eti[u] + 




g'^dudt + 


f‘‘r^{l + r‘^)drdu]dt 


^[Ti,T2\ 


) 


Proof. We have by (21) that 

Div ('^J) = {nAdSU + 2m) dtu = fdtu. 

Integrating this over S[Ti,T 2 ] applying the divergence theorem we pick up terms on the 
left hand side from , T 1 T 2 and . A straightforward calculation shows 
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We also find 


^ = [ r^{l + r^) (dtu) {dru)dujdt 


r"‘{dtuy + r‘‘{l + r‘^y{dru) 


^dt{ru) + r(l + y){dru)'^ 


dojdt. 


As r —?• oo, the term in braces may be replaced with g, so we have 


lim / 

1 — >0O Jj^lTi,T2] 


1 




2 hlTl,T2] 


^{dtuy + ^{druy — g^j dujdt 


Applying Lemma 2.1, we have: 
Et [ttl — Et. [u] = - [ 


—r‘^{dtuy — r^{druy + g^) dcodt + 


fdtur^drdtduj 


’[Tl.t] 


SO that 


(61) 

Now 


sup EAu] < EtAu] -\— [ g'^dujdt+-— -- [ ^ * \ ydrdtduj 

+ {T 2 —Ti) j + r‘^)drdu}dt. 


1 f jdtuy 

4(r2-ri)4^^ ^^j i + r2 


ydrdtdu < - sup 


{dtuf 


t£[Ti,T2] JT.t '^+r‘ 


ydrduj < - sup EAu]. 
2 te[Ti,r2] 


Applying this estimate to (61) and absorbing the energy term on the left hand side, we are 
done. □ 

After taking A to be large, we now construct true solutions of the wave equation with 
dissipative boundary conditions which only have a small loss of energy. To do this, we 
define ue,l,\ to be ue,l\ (defined by (55) in the previous subsection) multiplied by a 
constant factor in such a way that 

Eo{ue,l,\) = 1- 

Let Ue,l,x solve the homogeneous initial-boundary value problem: 

^9AdS^E,L,X + ‘^Ue,L,X = 0 

subject to dissipative boundary conditions 

d{rUE,L,x) , 2^(^^e,l,a) 


dt 


+ r 


dr 


0 , 


as r —>■ oo, 
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and initial conditions 

UE,L,\\t=Q = UE,L,\\t=Q^ 
We obtain the following theorem: 


dUE,L,\ 

dt 


t=o 


duE,L,X 


dt 


t=0 


Theorem 5.7. Fix E > 0 and 0 < |L| < E. Fix also T,e > 0. There exists a solution 
Ue,l,x of the homogeneous conformally coupled wave equation with dissipative boundary 
conditions such that Ue,l,x has energy 1 at time 0, and 


inf Et{UE,L,x)> 
te[o,T] 


e-Ve^-l'^ \ 

E + VE^ - L2 j 


where C > 0 is some universal constant. 


Proof Consider Ue,l,x as defined above, which is a solution to the homogeneous conformal 
wave equation. By Lemma [5. 12 as A —)■ oo, we have 


and 


/ {Ug^^^UE,L,x + 2m£;,l,a)^t^(1 + r'^) dr dudt^O 
■^'S’[0,T] 

{dt{ruE,L,x) + r^dr{ruE,L,x)f duj dt 0. 


qo,T] 


Therefore, after taking A to be sufficiently large and applying Theorem 5.6 to Ue,l,x — 
ue,l,x-, we can assume 

sup Et[UE,L,x - ue,l,x] < e - 

i6[0,T] 

On the other hand, by the construction of ue,l,x, we have 

'E - VE^ - 


' , _\ 2Af+l 

E - y/E^ - L2 ^ ^ 

E + ^/W^^ J ~ 

for some C > 0. The results follow straightforwardly. 


inf EtluE,L,x] > 
ie[o,T] 


E + VE2 - L2 


□ 


In particular, by taking \L\ sufficiently close to E and e sufficiently small, this implies 
that on the time interval [0, T], the loss of energy can be arbitrarily small. This also implies 
that any uniform integrated decay estimates without loss do not hold: 


Corollary 5.8. There exists no constant C > 0, such that 

Et[u\dt < CEo[u] 

holds for every solution u to the conformal wave equation with finite initial energy subject 
to dissipative boundary conditions. 

Similarly, there exists no continuous positive function f : M+ —)• M+, such that /(t) —)• 0 
as t ^ oo and 

Et[u] < f{t)Eo[u], 
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holds for every solution u to the conformal wave equation with finite initial energy subject 
to dissipative boundary conditions. 


6. Generalizations 

6.1. Alternative boundary conditions. 


6.1.1. Conformal Wave. We assumed that our solution u satisfied the boundary condition 

0 , 


d{ru) ^ 2 ^(™) 


dt dr 

We can also consider the boundary conditions 

d{ru) 29{ru) 


as r —oo . 


+ (d{uj)r^ 


0 , 


as r —>■ oo . 


where /3 : 


dt dr 

is a smooth function satisfying the uniform positivity condition: 

/3{uj) > 


for some k > 0 for all w G S^. This can be treated exactly as above, with the same 
results although the constants in the various estimates will now depend on /3. One could 
also imagine adding some small multiples of tangential derivatives of u to the boundary 
condition. This can also be handled by the methods above, but this will require combining 
the energy and integrated decay estimates. 


6.1.2. Maxwell. The boundary conditions that we assumed, 

{Ea + eA^Hs) —)• 0, as r —oo, 

can also be generalised without materially affecting the results. In particular, we could 
choose as boundary conditions 

{Ea + PAB{oj)e^^He) 0, as r oo 

where the symmetric matrix valued function /? : —)• M(2 x 2) satisfies a uniform positivity 

bound: 


(62) 

for some k > 0 and for all ^ 
Leontovic boundary condition 


G 


a; G S^. In particular, our results hold for any 


87] satisfying (62). Again, one could also permit other 
components of E, H to appear in the boundary condition with small coefficients and this 
can be handled by combining the energy and integrated decay estimates. 
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6.1.3. Bianchi. Recall that we are considering boundary conditions of the form: 

(63) (^Eab + ^a^Hb)c) -^0, as r oo. 

To generalise these, let us introduce a 4—tensor on S^, (3abcd{oj), symmetric on its first 
and last pairs of indices and also under interchange of the first and last pair of indices: 


Pabcd = /3{ab)cd = /3ab{cd) = I^cdab 

we also require that /? is trace free on its first (or last) indices, i.e. /3^abc = 0. In other 
words, /3 represents a symmetric bilinear form on the space of symmetric trace-free tensors 
at each point of . We can consider boundary conditions 


(64) [pAB^^EcD + e{A^HB)c) -^0, as r oo. 

Provided that I3abcd is uniformly close to the canonical inner product on symmetric trace- 
free tensors, our methods will apply. More concretely, there is a 5 > 0, which could be 
explicitly calculated, such that if 


(1 -I- 5) > /3abcd{^)'^^^'^^^ > (1 — <5) 


holds for any symmetric traceless H and w G S^, then our results hold for the boundary 
conditions (64). 

Of course, this does not imply that boundary conditions which don’t satisfy this inequal¬ 
ity lead to growth, simply that our approach breaks down when the boundary conditions 
are too far from the “optimally dissipative” ones (63). 


6.2. The Dirichlet problem for (B). Let us briefly discuss the Dirichlet problem for 
the Weyl tensor and its relation to the Dirichlet problem from the point of view of metric 
perturbations. In the latter, one wishes to fix, to linear order, the conformal class of 
the metric at infinity. Fixing the conformal class to be that of the unperturbed anti-de 
Sitter spacetime is equivalent to requiring that the Cotton-York tensor of the perturbed 
boundary metric vanishes. One may verify that the Cotton-York tensor of the boundary 
metric vanishes if and only if 


r^HAB 




0 


as r —)■ oo , 


as can be established by considering the metric in Fefferman-Graham coordinates. Let us 
illustrate in what way fixing the conformal class of the metric on the boundary is a more 


restrictive condition than fixing Dirichlet-conditions on the Weyl tensor, 


r^HAB 


0 . 


It is possible to extract from the Bianchi equations a symmetric hyperbolic system on 
involving only v^Eat, where r^HAB appears as a sonrce term. Using this system 

it is possible to show that if t’^Eat and r^H^rr vanish at infinity for the initial data then 
this condition propagates. Moreover, it is easy to see how to constrnct a large class initial 
data satisfying this vanishing condition at infinity, as well as a large class not satisfying 
it illustrating that fixing the conformal class of the metric on the boundary is a more 
restrictive condition than purely fixing Dirichlet-conditions on the Weyl tensor. 
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In conclusion, for initial data satisfying the vanishing condition, we may return to the 
estimate (42) and establish directly that solutions of the Bianchi system representing a 
linearised gravitational perturbation fixing the conformal class of the boundary metric 
are bounded. This is in accordance with the results of [29j, in which it is shown that 
metric perturbations obeying the linearised Einstein equations can be decomposed into 
components which separately obey wave equations admitting a conserved energy. 


6.3. The relation to the Teukolsky equations. We finally contrast our result with an 
alternative approach to study the spin 2 eqnations on AdS, which has a large tradition 
in the asymptotically flat context and relies on certain curvature components satisfying 
decoupled wave equations. As we will see below, however, in the AdS context this approach 
merely obscures the geometric nature of the problem and does not provide any obvious 
simplihcation as the resulting decoupled equations couple via the boundary conditions 
(and moreover do not admit a conserved energy). 

To decouple the spin 2 equations, we take the standard 9, 4> coordinates for the spherical 
direction^^ and choose as basis ei = r~^d 0 , e 2 = (r sin0)“^cl0. We then write: 

= Ell ~ E22 T 2 i?i 2 T * (^^11 “ H22 i 2E12) ■ 


These quantities obey the Teukolsky equations: 

„2 1 + 

1 . / . 1 


|2iTi± 


0 = -- 

1 +r2 * 


±4- 

l+r2 * 


dr 


1 + r' 


,dr [^(l + r^)'I'^ 


(65) 


+ {sin 680 ^^) + 

sm 9 


sm 


^ sm 9 


sin^ 9 


- 2 T 


Once have been found, the other components of W can be recovered by solving an 
elliptic syst em coupled to the symmetric hyperbolic system in the boundary discussed in 
Section 6 ^ It might appear that one can simply study the decoupled equations for 
separately. Unfortunately, in general, the correct boundary conditions couple the equations. 

Let us see what the appropriate boundary conditions to impose on are in order to 
fix the conformal class of the boundary metric. Clearly the vanishing of Hab at ^ is 
equivalent to the condition 


( 66 ) 


*(T+ -T-)| 


0 


as r 


00. 


This however only gives us one condition for two equations. For the other condition we use 
the fact t hat i n the context of the Dirichlet problem (fixing the conformal class discussed 
in Section 6.2) we know that r^EAri'f'^Hrf vanish on the boundary]^ Inserting this into 


22i 


The {r,t,0,cf>) do not quite cover AdS, and so some care must be taken at the axis. For the purposes 
of this section, we shall ignore this difficulty. 

In the general case, their trace on the boundary can de determined by solving transport equations as 


23 


outlined in Section 


6.2 
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the Bianchi equations we can derive that 


.2 d 




0 as r —)> oo. 


This gives us a Neumann condition for the Teukolsky equations: 


(67) 


d 


rr^(T+ + T )1 

or 


0 as r —oo. 


We should of course not be surprised that the two Teukolsky equations couple at the 
boundary. The two scalar functions represent the outgoing and ingoing radiative 
degrees of freedom. Since the Dirichlet boundary conditions are reflecting, one should 
of course expect that the two components couple at the boundary. The pair of equations 


(65) with the boundary conditions (66), (67) forms a well posed system of wave equations. 

In any case, there seems to be no 


as may be seen for instance with the methods of 
advantage in studying this coupled system of wave equations over the first order techniques 
of this paper. 

We remark that the Teukolsky approach was recently used in |19) to consider pertur¬ 
bations of the Kerr-AdS family of metrics (which includes anti-de Sitter for a = m = 0). 
In this paper, the Teukolsky equation is separated and a boundary condition (preserving 
the conformal class of the metric at infinity) is proposed for the radial part of each mode 
of separately. This appears to contradict our discussion above. When one examines 
equations (3.9-12) of [19] one sees spectral parameters appearing up to fourth ordeip^ 
Returning to a physical space picture, these will appear as fourth order operators on the 
boundary. Accordingly, it is far from clear whether these boundary conditions can be 
meaningfully interpreted as giving boundary conditions for a dynamical evolution prob¬ 
lem. Indeed, our heuristic argument for the coupling strongly suggests that the conditions 
of [T9| understood as boundary conditions for a dynamical problem cannot give rise to a 
well posed evolution. That is not to say that these boundary conditions are not suitable 
for calculating (quasi)normal modes: any such mode will certainly obey these conditions, 
providing a useful trick to simplify such computations. 


7. Appendix 

7.1. Proof of Lemma 15.31 

Proof. Let us first consider the first two terms of K. We set 

Ki := dr (j'^dr (^u\/ 1 -|- (^^u) ca — {r^ A'^)dr Cr . 


94 

In the presence of rotation things are even worse, as the square root in the boundary conditions implies 
that the operator on the boundary is non-local. 
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Calculating with the expression for the divergence of a vector field (18), we find 

1 d 

Div iTi = dr (^r^dr (^u\/l + - r'^y/l + r‘^{r'^Au)dr (r'^^u 

= dr (r^dr (u\/l + — \/l + r'^ dr [r'^u) 


+ A dr (r^dr (uy/l + ^ ^ (r^^i 


dr . 


= dr (r^dr (u^/l + — -x/l + dr [r'^u) 


+ 


dr {r'^dr {r'^j^u\/l + 


(y'^u) Vl + r^ , ^ s d 
—- 1 : 2 — i^yAU) — 


dr - 


dr ir'^^ 


= dr (r^dr (u\/l + a^^u — + r^ dr (r'^u) 


+ 


V^TT^ 


(r^^u) dr (r'^^u) . 


Finally, taking 


we calculate 


K2 :=- 


2(1 +r2) 


Iryul^ Cr , 


Div K 2 = — 


]^d^ 

ry* 2 

r 

\/l + r" 


\/l + r2 




l 


2(1 + r2) 


r^u\ 


[r^P au) dr yr'P^uj — 
{r'Pj^v) dr [r'P^u^ — 


1 _|_ J.2 

2r2 

3 \ r'p 


^ \2 d 


dr y(l + r2)i 


u\ 


2(1 + r 2)2 

Adding these two contributions, we arrive at the result. 


□ 


7.2. Proof of Lemma 15.51 Let us introduce 

2 


a := — 


V^' 


dr ^r-\/r+^^^A^ 


Hr. 


From the expression (18) for the divergence of a vector field, we can quickly establish that, 
owing to cancellation between the mixed partial derivatives, we have 

Div K = 2\/l + r2^ (- 2^ (ry^l + rmB) '’ f 


dr V\/l + r2 


dr 


ViV' 


^9^ 

r2 dr Y^l + r2 


m 
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SO that 


a 




= Div K 2 


= Div K + 2 


dr VVl + 1-2 

1 + (9 / r'^Hr \ d , 2 tt \ ^ d 

^ V rj + 


dr \ -|- r2 




7-2 Y J dr 

2 


+ 


= 2 


2r 


Vl + 

Hr 

Vl + 

d 

^2 

dr 


\/l + r 


d 


d 




dr \ 


m 


(1 +r2)2 


\/l + r2 5r 
T 


+ 


(1 + r^) 2 


+ Div iiT. 


3 U-'f’ 


Here we have used the constraint equation to pass from the first to the second line. This 
completes the proof. 


7.3. Proof of Lemma 5.7\ This is a straightforward calculation with the formula for the 
divergence of a vector field (18). We hnd 


1 

Div i7 = ^ — (2rWl + r^f^HBcH^r - [r^HBr 
or \ 

f^r^^[r(l + rV®c] 


Vl + r2 dr 


2rHBr d 
\/l + r2 dr 
2rHBr d 


r2(l + + 2(1 + Hbc^ 

J cr [Vl +r2 

r^l + r^)1/^H^c\ - -^^L=^^BB-r^ [r(l + r‘^)ii^c\ 
J Vl + d'T 




Vl + r2 dr 
- 2 rHBr^ {r^H\) 


dr 

2(1 + r^)2 d 


2r2 


dr 


^HBr) 


^ ^ ttB 


Vl + r2 dr 


•H^r) - 




(1 + r2)2 


I®' I' ^ 


Here, we have used the constraint equation in passing from the second equality to the third 
by replacing Hbc with a term involving dr{r^HBr)- 
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